VOLUME AUGUST 1959 NUMBER 


Canadian 
Journal Physics 


Editor: DUCKWORTH 


Associate Editors: 


ELLIOTT, Atomic Energy Canada, Ltd., Chalk River 
FOSTER, McGill University 

HERZBERG, National Research Council Canada 

LEPRINCE-RINGUET, Ecole Polytechnique, Paris 

SARGENT, Queen’s University 

WATSON, University Toronto 

WOONTON, McGill University 


Published THE NATIONAL RESEARCH COUNCIL 


OTTAWA CANADA 


| 
: 


CANADIAN JOURNAL PHYSICS 


Under the authority the Chairman the Committee the Privy Council Scientific 
and Industrial Research, the National Research Council issues THE CANADIAN JOURNAL 
and five other journals devoted the publication, English French, the 
results original scientific research. Matters general policy concerning these journals are 
the responsibility joint Editorial Board consisting of: members representing the National 
Research Council Canada; the Editors the Journals; and members representing the 
Royal Society Canada and four other scientific societies. 


EDITORIAL BOARD 


Representatives the National Research Council 


McT. Cowan, University British Thode (Chairman), McMaster University 
Columbia Thomson, McGill University 
Gauthier, University Montreal 


Editors the Journals 


Bailey, University Toronto Elliott, Montreal Neurological Institute 
Cameron, Macdonald College Léo Marion, National Research Council 
Duckworth, McMaster University Murray, University Western Ontario 


Representatives Societies 


Bailey, University Toronto Elliott, Montreal Neurological Institute 
Royal Society Canada Canadian Physiological Society 

Cameron, Macdonald College Gendron, University Ottawa 
Royal Society Canada Chemical Institute Canada 

Royal Society Canada Canadian Society Microbiologists 
Canadian Association Physicists Thorvaldson, University Saskatchewan 


Royal Society Canada 


officio 


Léo Marion (Editor-in-Chief), National Research Council 
Marshall (Administration and Awards), National Research 


Manuscripts for publication should submitted Dr. Duckworth, Editor, Canadian 
Journal Physics, Hamilton College, McMaster University, Hamilton, Ontario. 
For instructions preparation copy, see NOTES CONTRIBUTORS (back cover). 


Proof, correspondence concerning proof, and orders for reprints should sent the Manager, 
Editorial Office (Research Journals), Division Administration and Awards, National 
Research Council, Ottawa Canada. 


Subscriptions, renewals, requests for single back numbers, and all remittances should sent 
Division Administration and Awards, National Research Council, Ottawa Canada. 
Remittances should made payable the Receiver General Canada, credit National 
Research Council. 


The journals published, frequency publication, and prices are: 


Canadian Journal Biochemistry and Physiology Monthly $9.00 year 
Canadian Journal Botany Bimonthly $6.00 year 
Canadian Journal Chemistry Monthly $12.00 year 
Canadian Journal Microbiology Bimonthly $6.00 year 
Canadian Journal Physics Monthly $9.00 year 
Canadian Journal Zoology Bimonthly $5.00 year 


The price regular single numbers all journals $2.00. 


| 
j 


Canadian Journal Physics 


Issued THE NATIONAL RESEARCH COUNCIL CANADA 


AUGUST 1959 NUMBER 


VOLUME 


COMPARISONS BETWEEN FREE ENERGIES ACTIVATION FOR 
GRAIN GROWTH, GRAIN-BOUNDARY SELF-DIFFUSION, AND 
LIQUID 


ABSTRACT 


comparison made between the free energies activation for boundary 
migration during normal grain growth, grain-boundary self-diffusion, and liquid 
self-diffusion. The equality these, within the accuracy experiment and 
analysis, supports the assumption identical elementary processes the three 
cases. detailed discussion given regarding the various constants used the 
equations from which the free energies activation the above-mentioned 
processes are derived. 


INTRODUCTION 

The near equality between the values obtained for the activation energies 
and free energies activation for grain growth zone-refined tin and lead 
and those for liquid self-diffusion has recently been noted Holmes and 
Winegard (1959), who suggested that the elemental process may the 
same both cases. comparison the experimental results for normal 
grain growth with the theory Cole, Feltham, and Gilliam (1954) indicated 
that over the temperature range investigated, the mechanism boundary 
migration during normal grain growth involves single-atom process. Aust and 
Rutter (1959) have made similar observations their work single-boundary 
migration zone-refined lead. Holmes and Winegard (1959), applying the 
theory the kinetics solidification given Jackson and Chalmers 
(1956), have further suggested that the energy barrier overcome 
transferring atom from the solid the liquid during melting comparable 
the activation energy for grain growth high-purity zone-refined materials. 
has also been suggested Beck (1950) and Turnbull (1951) that the 
elementary process for grain-boundary migration closely related the 
elementary process grain-boundary self-diffusion that the free energies 
activation for the two processes should equal. 

this paper, values the free energies activation for grain-boundary 
self-diffusion, liquid self-diffusion, and, where possible, boundary migration 
during normal grain growth for various metals have been derived for purposes 
comparison. 


received March 13, 1959. 
Contribution from the Department Metallurgical Engineering, University Toronto, 
Toronto, Ontario. 


Can. Phys. Vol. (1959) 


899 


2 
> 
| 


900 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


FREE ENERGIES ACTIVATION FOR LIQUID SELF-DIFFUSION AND 
GRAIN-BOUNDARY SELF-DIFFUSION 


Values the free energy activation for self-diffusion the liquid close 
the freezing temperature are shown Table The method Mackenzie 


TABLE 


Calculated free energies activation the melting temperature, using jump distance 


Liquid self-diffusion 


Calculated 


from 
Mackenzie Calculated Grain growth 
Grain-boundary and Hillig from zone- 
self-diffusion, diffusion experimental refined metals, 
Metal coefficients data kcal/g-atom 
Lead 6.6 4.8 6.2 
(Okkerse 1954) (Rothman and (Bolling and 
Hall 1956) Winegard 1958) 
Tin 3.4 
(Careri and (Holmes and 
Paoletti 1955) Winegard 1959) 
Silver 
Turnbull 1951) al. 1958) 
Zinc 6.8 
1954) 
Cadmium 4.7 


(Wadja al. 
1955) 


and Hillig (1958) was used obtain liquid self-diffusion coefficient the 
freezing temperature from which the appropriate free energy activation 
was calculated using the expression: 


where the jump distance, taken This assumes that self-diffusion 
the liquid near equilibrium process which transition state theory 
may properly applied. The form equation (1) also assumes the simplest 
general case six equivalent directions for the migration each atom 
discussed LeClaire (1949). The free energy activation for liquid self- 
diffusion was also calculated from experimental data where available, 
indicated Table 

Also listed Table are values the free energy activation for grain- 
boundary self-diffusion derived applying transition state theory available 
experimental data exactly the same manner previously indicated, and 
following Hoffman and Turnbull (1951) and Wadja al. (1955), assuming 

FREE ENERGIES ACTIVATION FOR NORMAL GRAIN GROWTH 
ZONE-REFINED METALS 

Turnbull (1951) has developed equation for grain-boundary migration 

using the formalism absolute reaction rate theory, viz. 
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where the rate grain-boundary migration; 
the frequency factor; 
interatomic distance the boundary; 
aconstant depending grain shape; 
specific grain-boundary free energy; 
the boundary radius curvature; 
the gram-atomic volume the metal; 
AS, the entropy activation per gram atom; 
and the measured activation energy per gram atom. 


The driving force for migration was assumed the reduction total 
surface free energy; and the other terms are treated constants, the 


equation may written 


From this, the equation grain growth, viz. 
(4) 


nature grain growth. proper choice constants equation (2), the 
entropy activation, and thus the free energy activation per gram atom, 
may obtained from 

(5) 


the work Bolling and Winegard (1958) grain growth zone-refined 
lead, the following values were chosen for the constants: 2.72, 
vation for grain growth zone-refined tin near the melting temperature was 
calculated Holmes and Winegard (1959) using the same constants, with 
the exception that for tin 16.2 and was taken 100 
estimated Aust and Chalmers. These calculated values free 
energy activation for grain growth zone-refined metals are also listed 


Table 
DISCUSSION 
evident from Table that equality, within the limits experiment 
and analysis, exists between the derived free energy activation for liquid 
self-diffusion near the freezing temperature and for grain-boundary self- 
diffusion and grain growth near the melting temperature. This might re- 
garded indicating that the elemental process may similar all three 


cases. 
felt, however, that critical discussion required concerning the 
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choice constants involved and the assumptions made deriving the free 
energies activation. 

The following discussion thus concerned with the constants used 
equations (1) and (2). 


(a) Jump Distance 

Following Burke and Turnbull (1952) somewhat arbitrary jump distance 
was chosen both for the boundary region and the liquid state for the 
values listed Table 

With regard the liquid state, there appear present two possi- 
bilities: 

(i) The jump distance equal the distance the first maximum the 
radial distribution function given X-ray measurements. this possible 
then would appear that the calculated entropy activation for liquid self- 
diffusion would have negative value many instances. This may mean 
that before jump could occur, the disordered liquid would have become 
ordered assembling its free volume into discrete hole molecular 
dimensions, suggested Nachtrieb (1957). Only then would single atom 
jump interatomic distance the vicinity its stationary neighbors. 

(ii) The entropy activation zero that the transition state theory 
applies the jump distance must less than considered above. this case, 
diffusion not accomplished isolated jumps single atoms, but the 
co-operative readjustment many atoms. These possibilities have been 
discussed Nachtrieb (1957). 

The arbitrary value chosen, viz. somewhat between the two extremes. 


Boundary Width 

The value derived for the free energy activation for grain-boundary 
self-diffusion depends very much the choice boundary width. This 
discussed some detail McLean (1957), who concludes that for ordinary 
high-angle boundaries choice three atomic diameters reasonable. 
this value chosen rather than then the calculated values the free 
energy activation for boundary self-diffusion would increased somewhat, 
seen Table II. 


(c) The Constant 

This constant depends grain-boundary shape and, according Burke 
and Turnbull (1952), should the order Feltham (1957) suggests 
that this constant should for cylinder, for sphere, and concludes that 
the value 3/2, whereas Bolling and Winegard (1958) chose, arbitrarily, 
the lowest value, unity. Using the factor 3/2 for and assuming the other 
constants before, the free energies activation for grain growth tin and 
lead near their melting temperatures are 6.9 kcal/g-atom and 3.7 kcal/g-atom, 
respectively, rather than 6.2 and 3.2 kcal/g-atom obtained when assumed 
unity. 


(d) The Specific Grain-Boundary Free Energy, ergs 
The values 200 ergs for lead and 100 ergs for tin, estimated 
Aust and Chalmers (1950, 1951), were chosen for calculating the values 
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TABLE 


Metal Jump distance used, 
Lead 5.70 —3.92 3.4 
3.25 1.272 
Tin 4.25 —2.5 
3.00 1.716 
Silver 0.78 
10.3 2.0 


(b) Grain growth close melting temperature 


Free energy Entropy Values constants, etc. used, 


Metal keal/g-atom cal/g-atom ‘a’ taken unity 


(c) Grain- 


and silver close melting temperature 


boundary self-diffusion lead 


Free energy activation, 


Metal kcal/g-atom Jump distance, Boundary width, 
Lead 5.56 1.272 5.0 

9.10 3.4 5.0 

6.60 2.0 5.0 

4.66 1.272 2.0 

5.80 3.4 2.0 

7.50 2.0 10.0 
Silver 8.1 0.78 5.0 

12.70 2.0 5.0 

14.1 2.0 8.7 


shown Table The assumption has been made that these values are 
constant over the temperature range investigated. Any variation with tem- 
perature these values will mean that the entropy activation for grain 
growth also function temperature. 

high-angle boundary may regarded supercooled layer few 
atoms thick, then the grain-boundary free energy may expressed the 
sum the excess free energy the supercooled layer over the solid, plus 


(a) Liquid self-diffusion close freezing temperature 
Free energy activation activation 
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twice the liquid-solid interfacial free energy. This model, discussed 
McLean (1957), requires assumption regarding the effective width the 
boundary. One might assume, however, that near the melting temperature 
the values and 109 ergs viz. twice the interfacial free 
energy given Turnbull (1950) for lead and tin, respectively, might repre- 
sent the specific boundary free energy. The consequences this assumption 
are made evident Table II. 


(e) The Constant 

The formal theory recrystallization kinetics, discussed Burke and 
Turnbull (1952), describes the process isothermal recrystallization 
terms the parameters the nucleation frequency, and the linear rate 
growth. Burke and Turnbull (1952) point out that there good evidence 
that the kinetics normal grain growth after completion recrystallization 
are controlled large part parameter that closely related the 
values measured primary recrystallization. account for the fact that the 
normal grain growth cannot have the same formal relation atomic 
mobility recrystallization, the geometric factor ‘a’ introduced. This 
factor discussed Burke and Turnbull (1952), who conclude that 
general its value less than unity. The value used the calculations the 
free energies activation this paper was arbitrarily chosen unity. 

view the above considerations, interesting see just how the 
choice jump distance, boundary width, and choice constants used affect 
the calculated values the free energies activation. 

Table shows possible values for the free energies and entropies 
activation for liquid self-diffusion lead, tin, and silver. Table lists both 
the free energies and entropies activation for grain growth zone-refined 
lead and tin close their respective melting temperatures, and Table 
represents grain-boundary self-diffusion lead, using the experimental data 
Okkerse (1954), and silver, using the results Hoffman and Turnbull 
(1951). 

From the calculations Table evident that, the boundary 
region regarded being definite region width and that within this 
region there exist six equivalent directions for the migration each atom, 
then near the melting temperature near equality exists between the calcu- 
free energies activation for grain growth, grain-boundary self-diffusion, 
and liquid self-diffusion, assume that similar jump distance applies 


all cases. 


TEMPERATURE DEPENDENCE 
might argued that the comparisons made Tables and are valid 
only near the melting the case lead, tin, and silver, 
possible make comparisons over range temperature. 
Table indicates that the temperature variations the free energy 
activation with regard grain growth and self-diffusion extrapolated 


below the freezing temperature both lead and tin are similar nature. 
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TABLE III 


grain-boundary for liquid for 
Temperature, self-diffusion, self-diffusion, grain growth, 
Lead 600 6.6 4.4 6.2 
550 4.4 6.2 
500 4.5 
400 8.6 4.5 6.3 
Tin 500 3.3 
460 3.4 3.4 
440 3.4 3.6 
Silver 1233 12.7 10.3 
1000 14.4 9.9 
15.1 9.8 


900 


However, would appear that the free energies activation for grain- 
boundary self-diffusion increase more rapidly the temperature decreased 
than found the case for both grain growth and liquid self-diffusion. 
However, the measurements grain-boundary self-diffusion were undertaken 
relatively impure materials and this difference the temperature depen- 
dence the free energy activation for grain-boundary self-diffusion may 
due increasing effect the impurities the boundary the tempera- 
ture lowered. critical assessment impurity effects must delayed 
until the conclusion experiments now progress concerning the effects 
specific impurities grain growth zone-refined metals. 


SUMMARY AND CONCLUSIONS 


Investigations normal grain growth zone-refined lead and tin have 
indicated equality between the values derived for the free energies 
activation over range temperatures for grain-boundary migration and 
liquid self-diffusion extrapolated below the temperature. 

Values the free energy activation for grain-boundary self-diffusion 
have been calculated for several metals, and near equality has been found 
between these values close the melting temperature and those calculated 
for liquid self-diffusion provided that the various assumptions made are not 
unreasonable. difference the variation these free energies activation 
with temperature was noted which may due the fact that the grain- 
boundary diffusion data was obtained using relatively impure materials. 
hoped that this will stimulate further work self-diffusion liquid metals 
refined materials. Despite the fact that the problem the proper choice 
constants, jump distances, grain-boundary thickness, difficult one, 
the above results suggest that the activated state for the process whereby 
atoms move liquid self-diffusion the same the case movement 


Comparison between free energies activation AF, over range temperature, using 
for 
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atoms within the grain-boundary region. This suggestion appears reasonable 
view our present state knowledge regarding the structure high-angle 
boundaries. 
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EFFECTIVE AND THERMAL NEUTRON CAPTURE CROSS 
SECTION AND RESONANCE CAPTURE INTEGRAL 


ABSTRACT 


The effective neutron capture cross section, the 2200 meters/sec neutron 
capture cross section, and the resonance capture integral, Pr-143 have 
been measured activation method. The values are: 93+10, 
and 190+25 barns respectively. The neutron fluxes have 
been measured with cobalt monitors. All these values have been calculated 
according the Westcott convention, which has been outlined insofar can 
applied activation work. 


INTRODUCTION 


order obtain quantitative understanding the changes reactivity 
which occur during the irradiation fuel nuclear reactor, knowledge 
the neutron capture cross sections the fission products necessary. 
present, the cross sections the more important fission products are known, 
but there are some fission products which may significant poisons for which 
long-lived, 13.6-day half life, and formed high yield fission, 
fission. Furthermore, knowledge the neutron capture cross sections 
the fission products necessary for precise determination fission product 
therefore interest measure the neutron capture cross section 
Pri, 

Preliminary results this work, reported recent review paper 
cross sections certain fission products Eastwood a/. (1958), lead 
value 124 barns for the effective cross section Since then the two 
components the effective capture cross section have been measured: that 
due the Maxwellian component the neutron flux, i.e. the thermal cross 
section, and that due the epithermal component, the so-called resonance 
capture integral. Also the experiments reported here have led effective 
cross section which lower than that reported earlier about barns; this 
change the preliminary value due reassessment the efficiency the 
counter used the earlier experiments. 


EXPERIMENTAL 
(a) Description the Activation Method Used the Measurement 
The capture cross sections were determined activation method 
involving the measurement the 17.3-minute formed the 
reaction The nuclear reactions leading the formation and 
upon irradiation are given Fig. seen that Pr' can 
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Cel42 (11 OF NATURAL Ce) 


Cel43 Pri43. ———» Nd143 
HOURS 136 DAYS 
o (ny) 
> Pri44 (STABLE) 
290 DAYS MIN 


Fic. nuclides mass 143 and 144 neutron irradiation Ce-142. 


formed two separate routes. relatively short irradiation large amount 
cerium oxide was made, and the formed was allowed decay out 
the reactor for days. Separate experiments, made according the 
methods described Sections and (c), showed that was present 
the material that time. Thus the route 


was unimportant for the formation This result was not unexpected, 
since the neutron cross section only barns (Roy and Yaffe 1956). 

Following the out-of-reactor period (during which more than 99% the 
formed the first irradiation had decayed aliquots the 
mixture were reirradiated for minutes, and the product examined 
quickly for This two-stage irradiation, without intervening chemical 
separation, clearly allows quantitative distinction made between the 
two possible routes for the formation from Formation 
must come from neutron capture and not from the alternate route 
involving neutron capture 

these experiments the measurement the activity rests the 
fact that the 3.0-Mev beta rays can counted the presence the 
intense activity absorbing the 0.92-Mev beta rays the latter 
aluminum. The amount formed was obtained from its beta disinte- 
gration rate and, after the had decayed, the amount used the 
experiment was also obtained from its beta disintegration rate. This method 
has the advantage that the amounts and Pr' were obtained from the 
same source that the determination chemical yield not necessary. 


Irradiation and Chemical Purification 

Ten aliquots the mixture were irradiated empty fuel rod 
position the NRX reactor. The irradiations were done groups two; 
for each group, the first sample was irradiated iron capsule; the second 
sample was irradiated iron capsule lined with 0.076 cadmium. For 
each group irradiations the power the reactor was constant. Immediately 
following each irradiation, cobalt wires, 0.0127 diameter, length, 
and 1.151 weight (Jervis 1957), were irradiated iron and cadmium- 
lined iron capsules from minutes order measure both the thermal 


and epithermal fluxes. 
After the short irradiations fast separation from was made the 


a 
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iodate method. The reirradiated cerium oxide was dissolved boiling con- 
centrated sulphuric acid which 30% hydrogen peroxide was added. 
After addition carrier, both and were precipitated with ammonia 
and the precipitate was dissolved nitric acid. The cerium was oxidized 
the tetravalent state the addition solid potassium bromate and was 
precipitated the iodate. The precipitate was centrifuged and discarded. 
The supernatant, containing the Pr, was neutralized with NaOH. the 
first few experiments the praseodymium precipitate was dissolved and two 
more cerium scavengings were made, but was found that the additional 
cerium scavengings did not improve appreciably the radiochemical purity 
precipitate. Therefore the subsequent experiments, these extra 
scavengings were omitted save time. The final precipitate was 
deposited onto filter paper. The filter paper was transferred 
aluminum disk and dried. The was now ready for counting. 


(c) Counting the Samples and Absolute Rate Determinations 

The disk which the precipitate was mounted was covered with 
400 mg/cm? aluminum. This thickness will absorb all the 0.92-Mev beta 
rays The source was counted with end window methane flow 
counter. The activities the sources were measured until the had decayed 
completely and also several hours later estimate the contribution 
long-lived activity. The efficiency the counter under that arrangement was 
determined counting aliquots source separated from the fission 
product beta proportional counter and the flow counter under 
the conditions just described. was found 5.8%. 

From the efficiency the counter, and the time elapsed between the end 
the irradiation and the beginning the count, the absolute disintegration 
rate was obtained. half life 17.3 minutes for was assumed 
(Peppard al. 1957). 

After the 17.3-minute had decayed, the was dissolved nitric 
acid and the solution was made known volume. Aliquots the solution 
were evaporated gold-coated films about for absolute 
beta disintegration rate beta proportional counter. was found 
necessary cover the sources with additional VYNS film coated with 
gold get sufficient conductivity avoid losses the counting rate. Two 
sources were made from each solution and their decay was followed for 
two half lives. The number atoms the reaction producing 
was obtained from that determination extrapolating the disintegration 
the midtime the irradiation using 13.6 days for the half life 
(Peppard 1957). 


(d) Identification 

all the determinations, the short-lived activities after subtraction 
residual activity had half lives ranging from 17.3 18.9 minutes. Figure 
shows the decay curves the activities for the group experiments labelled 
Table separate experiment was done show that the short-lived 
activity was due and not some other radioactive species carried 
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MINUTE 


(COUNTS PER 


ACTIVITY 


50 100 150 200 250 
TIME (MINUTES 
Fic. curves activities for the group experiments labelled Table IT; 
irradiation done without cadmium filter; irradiation done with cadmium filter. gross 
activity; activity due Pr-144; residual activity. 


down with the precipitate. sample cerium oxide, free from 
was irradiated for hour high-flux position the reactor. such 
irradiation the amount formed much too low lead detectable 
amount reaction. Praseodymium was separated and counted 
described Sections (6) and (c). activities with half lives the range 
minutes few hours were detected. 

The residual activity the precipitate found after the decay 
the reaction; has 2.1-Mev beta ray and would have 
been partly produced the decay and could have been present 
initially impurity the cerium oxide. 


5 ° 
10 
° 
° 
° 
° 
° 
O 
° 
105 
e 
q 
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(e) Flux Measurements 

The gamma activities cobalt were measured placing the monitors inside 
well-shaped hole crystal coupled RCA 5819 photo- 
multiplier tube. The pulses above 0.66 Mev were registered single-channel 
analyzer. The spectrometer was calibrated with Co® source whose absolute 
activities were determined and y-coincidence counting. The measured 
y-activities have been corrected for the self-shielding effects arising from the 
absorption and the scattering the neutrons the monitor. For the monitor 
used the experiments, Eastwood (1959) has found that the observed re- 
action rates without and with cadmium filters have multiplied 
1.04 and 1.62 respectively. These corrections will give the reaction rates that 
one normally would observe using infinitely thin cobalt wire. For the 
irradiations without cadmium filters, absorption and scattering the neutrons 
are responsible for the correction about the same amount; for those under 
cadmium filters, the large correction arises from the scattering resonance 
for neutrons 132 ev. The disintegration rates extrapolated the end 
the irradiation were obtained using half life 5.28 years for (Hughes 
and Schwartz 1958). 

The precision the measurements the thermal and epithermal fluxes 
was and respectively. 

(f) Errors 

The combined uncertainties obtaining the absolute disintegration rates 
and are estimated about 10% for the irradiations without 
cadmium and about for those done under cadmium filter. they are the 
major source errors, the final results are estimated precise within these 
limits, that 10% for the capture cross sections and 15% for the resonance 
capture integral. 

RESULTS 

obtain the flux from the cobalt monitor measurements and the capture 
cross sections the convention described Westcott (1958) was 
used. Eastwood (1958) have already described and applied this convention 
the measurement neutron fluxes and cross sections activation 
method. view the novelty the convention, however, may useful 
repeat its main features and give the equations applicable the results 
found activation work. 

This convention based equating the measured reaction rate, 
per atom capturing nuclide the product the effective cross section, 
and the reactor thermal neutron flux, where the total neutron density 
and 2200 meters/sec. The term usually called 
conventional flux. Thus 


The units are atoms product per unit time per atom capturing 
nuclide. For convenience activation experiments can given 


(2) 


\ 
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where the disintegration rate the product the end the irradia- 
tion, the number atoms capturing nuclide, the disintegration 
constant the product, and the time irradiation. 

well-moderated thermal neutron reactor, such NRX, the neutron 
spectrum consist two components: the first Maxwellian 
distribution with characteristic temperature without energy limits, 
the second flux distribution (the epithermal region) subject 
suitable lower energy limit but overlapping the Maxwellian component. 
The value thought about 0.125 for heavy water moderators. 
activation measurements, however, the distribution cut off 
0.5 the use cadmium filter. therefore more convenient use 
the effective cadmium cutoff rather than the lower energy limit. This 
interpretation makes the resonance capture integral smaller for 
absorber which follows the law the 0.125 0.5 region. 

Subject the assumptions Section the effective cross section 
defined the relation: 


(3) 


component the neutron spectrum, and and are factors which depend 
the departure the cross section from the law the Maxwellian 
and components respectively. 

The value characteristic the reactor and addition varies 
within the reactor. related the cadmium ratio, equation 
given later. When the neutron spectrum pure Maxwellian The 
values and have been tabulated for the most important nuclides, over 
range neutron temperatures (Westcott 1958). However, these tables cannot 
used activation work because the neutron temperatures are not usually 
known. The additional assumption that has made; fortunately 
this the case for most substances including the important detector 
obtain without measurement the neutron temperature, new factor, 
Sy, introduced such that 


(4) 


The value can obtained from equation 
report CRRP-785 (1958) which is: 


where the resonance integral the absorber, including the part, 
0.025 ev, and the cadmium cutoff energy taken 
0.5 for 0.076 cadmium. 

Equation (3) can now written 


/ 
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Following the convention, related the cadmium ratio, Reg, 
and equation Westcott a/. (1958), which for convenience can 
given the form: 

The cadmium ratio the ratio the reaction rate the absorber without 
cadmium the reaction rate with cadmium filter. The factor coefficient 
calculated from the variation the cadmium cross section with neutron 
energy and the thickness cadmium used the irradiation. Values 
for thicknesses cadmium common use are tabulated Westcott al. 
(1958); for 0.076 cadmium, the thickness used these experiments, 
with isotropic incidence neutrons, 2.2. 

From equations and reference standard (cobalt these 
experiments) and the resonance capture integral for any absorber can 
measured activation method. 


(a) Westcott’s Convention Flux Calculations 

the application the convention the calculation neutron fluxes 
from cobalt monitors the factor evaluated from equation From the 
measured cadmium ratio cobalt for the irradiation position used and the 
can deduced for cobalt using equation From this nvo follows from 
equation the reaction rate without cadmium filter obtained from 
equation 

The conventional neutron fluxes obtained this way and the data necessary 
for their calculation have been assembled Table The value for 


TABLE 


Conventional neutron fluxes measured with cobalt monitors 


Dura- 
tion of 
Experiment irradia- 
- ——— tion, dN/dt,*dN/dt,t R§ Rca. nvo, 
Group No. seconds d.p.s. d.p.s. N.3 sec"! Co-59 (T/To) barns n cm?-sec 


vo C 2 3,3¢ 3.9% 18 2.5 9 
A 1 NoCd 120 13,390 13,930 1.176X10 2.3710 19.4 1.80 0.0238 38.2 620x100 


120 442 717 10-10 


*Disintegration rate of Co® at the end of the irradiation, uncorrected for self-shielding. 
tDisintegration rate of Co® at the end of the irradiation, corrected for self-shielding. 
tNumber of atoms of Co present during the irradiation. 

§Rate of the Co (n,7)Co® reaction. 
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taken be: 36.68 barns, the average the two values BNL-325, 
2nd edition (Hughes and Schwartz 1958). The value the resonance capture 
integral including the 1/v part, for infinitely thin wire cobalt 
75+5 barns (Johnston 1959; Eastwood 1959). 


(b) Capture Cross Sections and Resonance Integral 

The effective cross section the reaction calculated 
from its observed reaction rate without cadmium filter and from the con- 
ventional flux using equation The flux has already been measured; the 
reaction rate calculated from equation 

The cadmium ratio for has been measured each group experiments. 
Then the factor for can calculated from equation where 
which constant the reactor for the irradiation position used, has already 
been determined. The 2200 meters/sec cross section can now calcu- 
lated from the values so, and application equation 

The neutrons the flux distribution have equal flux per energy 
interval well-moderated thermal neutron reactor such NRX 
that the resonance capture integral, can defined 
with the integral evaluated from Ecq 0.5 fission energies. Following 
the convention the resonance capture integral was obtained from equation 

The values and for and the data needed obtain them have 
been compiled Table The average the five determinations, including 


TABLE 


Capture cross section and resonance capture integral 


nvo 
Experiment Duration of (from 
irradiation, Table I), dN/dt,* Rea, a0, RI, 
Group No. seconds n/cm?-sec d.p.s. NT sec"! Pr-143 barns so barns — barns 
4 600 § 20% 10% 19,250 9.56%108 6 09% 
Qe 
B 3 1500 6§ 05 1038 37,420 10.03K102 5 34% 107% 
16.3 88.3 2.18 83.8 200 
( 5 1500 6.16% 10" 71,050 17.80X%10% 5.71 K10% 
185 92.7 
1500 107% 106% 50,550 12.25% § 107% 
20 .¢ 6 
} 4 1500 97 X 86,730 23 30K%10% 32K10 
| 80.1 2.13 84.8 
1500 1,545 20 63 x*10% 1IKIO 7 
Mean 03.0 80.0 190 
*)isintegration rate of Pris at the end of the irradiation 
TNumber of atoms of Pr present during the irradiation 
tRate of Pr(n.7) Pr reaction 


the estimated errors discussed Section Experimental, gives values 
Note that, unlike the observed the observed 
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Pr' have not been corrected for the self-shielding the neutrons the 
target material. assumed that self-shielding corrections are necessary 
for the amounts cerium oxide used the measurements, because the thermal 
absorption cross section cerium only 0.73 barn and the total cross section 
fairly constant barns from 1.0 1.0 Mev, with the exception 
small resonance 140 kev (Hughes and Schwartz 1958). 
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THE CUMULATIVE FISSION YIELDS LIGHT MASS 
FRAGMENTS THE THERMAL NEUTRON FISSION 


ABSTRACT 


The relative cumulative yields light mass fragments rubidium, 
strontium, yttrium, zirconium, molybdenum, and ruthenium formed the 
thermal neutron fission have been determined with mass spectrometer 
using the isotope dilution technique. 

The yields the fission product isotopes these elements, together with the 
previously published relative yields three krypton isotopes, have made possible 
the determination the absolute yields light fragment mass chains. 


INTRODUCTION 

Published values for the absolute cumulative yields the light mass frag- 
ments formed the thermal neutron fisson consist only few 
radiochemically determined vields summarized Katcoff (1958), the absolute 
vields and reported Krizhanskii and Murin (1958), and the yields 
the krypton isotopes determined the same time the present work 
Fritze al. (1958) using mass spectrometer. 

The absolute yields the strontium isotopes reported Krizhanskii and 
Murin (1958) are based the normalization their relative yields the 
5.98% yield Nd' reported Wiles al. (1956). These values should 
therefore renormalized the more accurate value 4.56% for the yield 
(Fickel and Tomlinson 1959). 

The cumulative yields mass chains the light mass fragments have 
been determined the present work and compared with the previously 
reported values. 


EXPERIMENTAL 

Plutonium was irradiated the NRX (Chalk River) reactor both 
sealed quartz tubes and alloy. The irradiation data are shown 
Table Complete dissolution the after irradiation was accomplished 
conversion the fluoride alternate treatment with and 
solutions during period several days. 

The relative yields the isotopes rubidium, strontium, and yttrium were 
determined using the mass spectrometer with prior chemical separations 
from the other fission products. However, carrier-free separation techniques 
were required the analyses molybdenum, and ruthenium 
fission products. 

the case zirconium, was taken its colloidal-forming 
tendency neutral separation was accomplished the dialysis 

Manuscript received May 19, 1959 

Contribution from Chemistry Department, University, Ontario 
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TABLE 
Irradiation data for plutonium samples 


Weight, Chemical 


Irradiation Average neutron flux 


Sample form time, days n/cm? sec 
0.01012 51.57 0.346 
0.0758 Pu-Al alloy 0.194 
7 0.10551 10°% Pu-Al alloy 30 6.6* 


*Only approximate, obtained from Reactor Operation Report. 


neutral solution the irradiated sample. described Schubert (1949), 
maximum efficiency was achieved greater than from solution 
completely free from such anions F~, and 

Excess acid was removed from the solution the irradiated sample 
repeated evaporations dryness and dissolutions water. The solution was 
dialyzed and the contents the dialyzer was analyzed with the mass spectro- 
meter. Whereas the separation from the other fission products was essentially 
complete this method, only partial separation from plutonium was 
obtained. 

The molybdenum was separated three extractions from aqueous 
solution using equal volumes ethyl ether saturated with 
The ether extracts were combined and analyzed with the mass spectrometer. 

Ruthenium was separated from the other elements glass apparatus 
used oxidize the ruthenium the volatile ruthenium tetroxide which was 
collected the distillate reduction with dilute solution and 

Samples and were used for obtaining absolute yield data. For samples 
and small fraction the solution was with the mass 
spectrometer for the presence contamination the fission product elements 
stable isotopes. Samples and were portions the same piece un- 
check was made for sample This sample was dissolved directly mixture 
the isotope diluents required tor the analysis. The method tsotope dilution 
has been described previously Petruska 


EXPERIMENTAL RESULTS 
The relative the isotopes the elements strontium, zirconium, 
yttrium, molybdenum, and ruthentum are shown Tables The data 
measured mass spectrometric ratios the strontium 


isotopes together with the ratios and the 


corrections for and have been made 
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Relative yields strontium isotopes the thermal neutron fission 


Time from end Ratio corrected: 


Mass 
Sample days Isotope ratio For contamination For decay 

0.644+0.013 
1.000 1.000 


*t4 for Sr8* = 51.8 days (Fickel and Tomlinson unpublished). 
tty for 27.7 years (Wiles al. 


Relative yields zirconium isotopes the thermal neutron fission Pu? 


Ratio corrected: 
irradiation, 
Sample days Isotope ratio For contamination 
3.361 
1.167+0.023 
1.000 1.000 
2.261+0.024 1.125+0.024 
0.1818 
93 1.000 1.000 
431 


*14 for Zr% 
tty for 


Time from end 


TABLE 


Mass 


239 


= 65.8 days (Fickel and Tomlinson unpublished data). 
= 58 days (Kahn et al. 1955). 


For decay 


0.791+0.021 
1.000 


1.312+0.018 


Relative yttrium isotopes the thermal neutron fission 


Sample 


Time from end Mass 
irradiation, spectrometric 
Isotope days ratio 


Note: ¢4 for Sr8* = 51.8 days (Fickel and Tomlinson unpublished). 
ty for Y! = 58 days (Kahn et al. 1955). 


Ratio 
corrected 
for decay 


1.00 


9 
| 
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TABLE 
Relative yields molybdenum isotopes the thermal neutron fission 


Time from end Mass Ratio corrected: 
irradiation, 
Sample days Isotope ratio For contamination For decay 
1.008 
0.797+0.014 
100 1.059 1.000 1.000 
0.3498 
0.994+0.017 0.794+0.017 0.794+0.017 
1.335+0.015 0.832+0.015 
100 1.204 1.000 1.000 


*Assuming 65.8-day half-life for Zr9* and 35-day half-life for Nb% (Engelkemeir et al. 1951; Cork et al. 1953). 


TABLE 


Ratio 


Time from end 
irradiation, spectrometric corrected 
Sample days ratio for decay 
668 101 1.000 
664.1 102 1.013+0.014 1.013+0.014 
104 
106 
104 1.000 
106 7730+0 


Note: ¢4 for Ru! = 39,7 days (Wright et al. 1957). 
ty for Ru! = 1.01 years average of 1.0 year (Schuman ef al. 1956) and 1,02 years (Merritt et al. 1957). 


Samples and were used for the determination the relative yields the 
zirconium isotopes shown Table III. The natural contamination 
and vields was obtained from the abundance should 
noted that the decay times used for the correction the measured yields 
and were not the same. The decay time for represents the time 
until the mass spectrometric analysis. However, the case the which 
grows through its 58.0-day precursor, the decay time only includes the 
time until the dialysis the sample. The relative yield was not deter- 
mined using sample since large decay correction requires accurate values 
for the half-life and also the time separation. Sample had decayed 
for over year; hence the time dialysis and half-life the precursor 
were not critical. 

The analysis the fission molybdenum, shown Table was complicated 
the small amounts natural molybdenum which were present the 
tungsten filaments. The measured relative abundances the fission molyb- 
denum isotopes were corrected for this natural molybdenum contamination 
from the amount mass indicated the mass spectrum. was also 


Relative fission yields ruthenium isotopes fission 


920 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


necessary correct the measured abundance for the amounts 
and which had not decayed the time analysis. For sample 
this correction was too large for accurate estimation the mass yield. 

One half sample together with sample was used the analysis 
fission ruthenium shown Table VI. The absence mass proved that 
natural ruthenium was present contaminant. The abundance was 
too small measured with any degree accuracy. The yield was 
determined relative using the isotope-diluted sample Since the 
estimation the relative yields from these measured abundances required the 
knowledge two half-lives, that not being well known, the final 
relative probably not accurate the mass spectrometric precision 
would indicate. 

addition the relative yields shown Tables VI, the relative 
the rubidium isotopes were determined. The measured ratio 
was 0.4517+0.0036. fission yield ratio 0.6354 0.0050 was obtained after 
correcting the for the decay the 10.27-year isomer during 
irradiation and the 0.3535 from the end irradiation the time 
analysis. was assumed that 23% the total mass chain decayed through 
the 10.27-year isomer while 77% decayed through the 4.36-hour isomer 
(Bergstr6m 1951). 

Table VII contains summary the relative yields the isotopes the 
various elements which were used the determination absolute fission 
vields. The values which were taken from Tables are averages weighted 
according the inverse square the standard deviations. 


TABLE VII 


Relative yields isotopes various light fragments the thermal neutron fission 


Relative yield Isotope Relative yield 

0.960 

020 
Zr” 1 O00 
1. 26640 022 


Samples and were utilized the determination the absolute yields. 
isotope dilution data along with the evaluated number atoms 
and per gram plutonium each sample are 


shown the relative yields shown Table VII 
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TABLE VIII 


Mass spectrometric and isotope dilution data for produced the thermal neutron 
fission 


No. atoms 


isotope Calculated 
Ratio before added per Ratio after fission yield, 
1.000 4.736 1.000 
2.202+0.017 1.828 2.068 
1.000 4.198 1.000 


TABLE 


Mass spectrometric and isotope dilution data for produced the thermal neutron 
fission 


No. atoms 


isotope Calculated 
Ratio before added per Ratio after fission yield, 
133 1.000 5.938 1.000 
133 1.000 7.853 
133 1.000 5.938 1.000 


137 


Mass spectrometric and isotope dilution data for produced the thermal neutron 


No. of atoms 


isotope Calculated 
Ratio before added per Ratio after fission vield, 
1 000 1 000 5 SHY 
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TABLE 


Mass spectrometric and isotope dilution data for produced the thermal neutron 
fission 


No. atoms 


isotope Calculated 
Ratio before added per Ratio after fission yield, 
isotope isotope atoms 10!8/g 
0.830+0.017 0.071 0066 13.50 
0.793+0.019 0.0541 0.5045+0.010 
0.830+0.017 0.079 15.60 


*Analysis of molybdenum solution used as isotope diluent: Mo” 0.932+0.002; Mo%’ 0.013+0.005; Mo** 
0.019+0.01. 


TABLE 


Mass spectrometric and isotope dilution data for produced the thermal neutron 
fission 


No. atoms 


isotope Calculated 
Ratio before added per Ratio after fission yield, 
Sample dilution dilution 
10) 1.000 13.60 


and the number atoms per gram plutonium for one isotope each fission 
element Tables the number atoms per gram pluto- 
nium has been tabulated columns and Table XIII for each the 


The relative the zirconium isotopes were normalized through the 
yield the mass chain whose value was obtained from the isotope dilution 
fission molybdenum. manner, the relative yields the krypton 
isotopes given were normalized the rubidium isotopes 


through the vield the mass chain which was determined isotope 


dilution 
For each the isotope-diluted samples, the number atoms pet 
plutonium was Since the absolute yield has been 


evaluated 6.90% and 1959), was possible normalize 


the atom vields the the other elements percentage 


The final values shown Table and graphically represented 


vields for masses and 112 1958) the 


vields from the interpolation and extrapolation plot are also 
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TABLE 


Cumulative fission yields the light fragments the thermal 
normalized the 6.90% yield 


Sample Sample Sample 
72-82 
0.29 
0.47 
1.460 0.5456 1.262 0.5251 0.535 
0.75 
2.487 0.9291 0.2150 0.8942 0.912 
1.438 3.819 1.440 3.164 1.368 1.43 
1.728 4.589 1.726 3.802 1.71 
0.3656 2.203 5.849 2.199 4.846 2.089 2.16 
2.59 
3.12 
3.94 
4.45 
14.97 5.630 12.97 5.61 
15.60 13.50 5.818 5.84 
6.44* 
1.882 7.020 7.05 
13.60 
13.78 5.94 
105 
107 
108 
109 
110 
112 
114 
115 
100 12 


yield 

*Interpolated values. 

tRadiochemical yields. 

DISCUSSION 

The relative vields obtained for the light mass fragments are less precise 
than those the previously relative vields the heavy mass 
fragments and 1959). Such the presence 
molybdenum tungsten filament corrections have 
precisions between have been obtained most cases. 

The accuracy the absolute vield values also depends the 
the relative the cases the rubidium, strontium, molybdenum, 
and this was isotope dilution. Phe 


the analyses shown Pable within tor most cases. zirconium 


924 CANADIAN JOURNAL PHYSICS. VOL. 37, 1959 


v 


YIELD 


FISSION 


a 6s 89 93 97 101 105 ws "3 
MASS 


Fic. mass yields thermal neutron fission 


was not isotope-diluted, but the relative yields were normalized those the 
other elements through the yield the mass chain obtained from molyb- 
denum. This procedure was used rather than direct isotope dilution because 
absorption characteristics this element solution make the equilibration 
the zirconium fission product and isotope diluent zirconium difficult. will 
noted from Table XIII that the yield 1.52 times greater than the 
yield The relative the mass and chains may also 
obtained from the relative yields the yttrium isotopes shown Table IV. 
The value obtained this manner but the precision only includes 
the uncertainty introduced from the mass spectrometric observations and not 
the systematic errors introduced the large half-life corrections. the factor 
1.44 had been used normalize the yield that the zirconium 
yields would have been about 5.5% lower than those given Table XIII. 
This inconsistency would substantially removed the yttrium ratios had 
been corrected for decay using instead for the half-life 

compare the yield obtained this work with the yield found 
Fritze al. (1958) necessary assume the branching ratio for the 
4.4-hour assumed that the short-lived krypton isomer 
decays directly rubidium, found that the 0.55% yield given Fritze 
al. becomes 0.539, which exact agreement with that found the present 
work. This has provided confirmation the branching ratio reported 
Bergstrém (1951), but slightly greater than the value estimated Blades 
and Thode (1955), who found that 78% the decayed directly 

Krizhanskii and Murin (1958) have published yields and formed 
the thermal neutron fission Both relative and absolute yields 
these two isotopes were significantly different from those found this work. 
These authors have given insufficient experimental detail assess the reason 
for the differences. 


8 
7 
6 
4 
‘4 
‘ 
‘ 
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The summation the yields the light mass fragments given Table XIII 
totals 100.13%, which assures confidence the values. The accuracy the 
absolute cumulative fission yields shown Table XIII considered 
better than 4%. 
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THE CUMULATIVE FISSION YIELDS HEAVY MASS 
NUCLIDES PRODUCED THE THERMAL FISSION 


ABSTRACT 


Mass spectrometric analysis cesium, barium, cerium, neodymium, and 
samarium produced the thermal fission has made possible the deter- 
mination the cumulative yields mass chains. With previously published 
mass spectrometric yields four xenon isotopes normalized the present data 
total cumulative yields are reported. 


INTRODUCTION 

Radiochemically determined yields nuclides formed the thermal 
neutron fission have been summarized Steinberg and Freedman 
(1951) and Steinberg and Glendenin (1955). The summation those for the 
heavy mass fragments assuming smooth curve totals only 93% indicating 
considerable error either the measured values those interpolated from the 
smooth curve. Mass spectrometric values for the absolute yields 
the heavy fragments were published Wiles and Tomlinson (1956) and later 
Krizhanskii a/. (1957). the work Wiles and Tomlinson (1956) only 
80% the fission product activity was recovered the procedure dissolving 
the irradiated was pointed out the discussion this work that the 
yields reported would inaccurate the losses the various elements were 
not equal. The yields obtained Krizhanskii a/. (1957), although essential 
agreement with those Wiles, can not considered completely independent 
and hence confirmation Wiles’ data, since only relative yields were actually 
measured and these yields normalized Wiles’ absolute yield values 
and 

The present work careful redetermination the yields reported 
Wiles al. together with further yield 


EXPERIMENTAL 

Two different samples were used determine the thermal neutron fission 
yields for plutonium. The first sample was 10.12 (96.1% 
sealed quartz capsule. This was irradiated along with flux monitor con- 
sisting 1.15-mg wire, 0.13 diameter, for 51.34 days flux 
thermal neutrons per square centimeter per second the NRX 
reactor, Chalk River. Eight months after the sample had been removed from 
the reactor was completely dissolved platinum crucible continuous 
hot treatment for approximately days. The quartz capsule was 
then washed several times with until all the activity had been removed. 
Without any chemical separation, small fraction (less than 1%) the solution 
was analyzed with mass spectrometer for the presence stable isotope con- 
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tamination each the fission product elements cesium, cerium, neodymium, 
and samarium. The remainder this solution was isotope-diluted using the 
method previously described Petruska (1955). The mass spectrometer 
was 10-in. radius, 90° sector, solid source instrument equipped with magnetic 
scanning and 10-stage electron multiplier assembly. 

second sample weighing was half alloy disk which had 
been irradiated for 23.88 days flux thermal neutrons per 
square centimeter per second the NRX reactor, Chalk River. The alloy 
contained 13.9% which 96.1% was order calculate the num- 
ber fissions this sample was necessary determine the neutron flux and 
the effective cross section. This was accomplished the following manner. 

The Pu-Al alloy was enclosed aluminum capsule about 2.5 
length. Cobalt wire flux monitors, 0.13 diameter and 1.56 and 1.66 
weight, were placed each end the capsule. The neutron flux was 
determined from the activity induced the cobalt which was measured 
well-type ionization chamber relative cobalt standards whose absolute 
disintegration rate had been previously determined absolute beta counting. 

The effective cross section the requires knowledge the neutron 
temperature. neutron temperature monitor, consisting samarium 
evenly deposited aluminum foil, was wrapped around the inside the 
capsule. 

The mass spectrometric ratio after the irradiation was used 
obtain the neutron temperature discussed below. 

Sixteen months after irradiation this sample was dissolved distilled 
The fission products were analyzed with the mass spectrometer both 
before and after isotope dilution the same manner for the sample. 
addition, analysis fission was carried out for this sample using 
solution enriched The concentrations the barium isotopes the 
enriched solution were determined means isotope dilution with 
known quantity natural barium. 


EXPERIMENTAL RESULTS 
The observed mass spectrometric ratios the fission products together with 
the relative fission yields the isotopes cesium, neodymium, cerium, and 
samarium are shown Tables IV, respectively, well previously 
reported values. Tables contain the isotope dilution data along with 


TABLE 
Relative yields the cesium isotopes produced the thermal neutron fission 


Relative fission yield 


Mass Literature 
spectrometric 
Isotope ratio This work Wiles (1956) Krizhanskii (1957) 
1.000 1.00 1.000 


*Corrected for 3.45 X 108-barn neutron absorption cross section of Xe! (Petruska et al. 1955). 
tCorrected for 26.6-year half-life of Cs? (Wiles et al. 1955). 
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TABLE 
Relative yields the neodymium isotopes the thermal neutron fission 


Relative fission yield 


Literature 
spectrometric 

Isotope ratio This work Wiles (1956) Krizhanskii (1957) 
1.000 1.000 1.00 1.00 

0.424+0.005 0.841+0.010 0.836+0.006 0.840+0.006 


*Correction for 334-barn neutron absorption cross section of Nd" was insignificant in this work (Walker and 


Thode 1953). 
tCorrected for 278-day half-life of Ce. 


TABLE III 
Relative yields cerium isotopes produced the thermal neutron fission 


Relative fission yield 


Mass Literature 
spectrometric 
Isotope ratios This work Wiles (1956) Krizhanskii (1957) 
1.12+0.01 1.12+0.01 1.11+0.02 0.85+0.01 
1.00 1.00 1.00 1.00 


*Corrected for 278-day half-life 
TABLE 
Relative yields samarium isotopes produced the thermal neutron fission 


Relative fission yield 


Mass Literature 

spectrometric 
Isotope ratios This work Wiles (1956) Krizhanskii (1957) 
1.00 1.209 1.209 
0.7410+0.0074 0.736+0.014 0.764+0.01 


*Corrected for 2.52-year half-life of Pm"? (Melaika et al. 1955). 
tCorrected for 12,800-barn cross section of Sm"™! (Melaika et al. 1955). 


TABLE 
Mass spectrometric isotope dilution data for produced the thermal neutron 


fission 


Isotope 
134 136 138 
Atoms added per gram 6.486 0.9485 1.940 
plutonium 
Relative abundances natural 0.242 0.781 7.166 
barium 
Mass spectrometric ratio after 1.000 
isotope dilution 
Atoms fission per gram 6.575 


— 
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TABLE 
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Mass spectrometric and isotope dilution data for produced the thermal neutron 


fission 


Isotope 
140 142 144 
Mass spectrometric ratios 
1.000 0.3694+0.0031 
1.000 
Ratio after dilution with natural 
cerium 
atoms per gram 1.000 00080 
plutonium 
atoms per gram 1.000 
plutonium 


Atoms fission per gram 
plutonium 


TABLE VII 


Mass spectrometric and isotope dilution data for produced the thermal neutron 


fission 


Isotope 
151 152 

Mass spectrometric ratios 

1.000 

1.000 0.864+0.010 
Ratio after dilution with natural samarium 

atoms per gram plutonium 1.000 13.95+0.40 

atoms per gram plutonium 1.000 
Atoms fission per gram plutonium 

1.638 


TABLE VIII 


*Corrected for 12,800-barn cross section of Sm! (Melaika et al. 


1955). 


Mass spectrometric and isotope dilution data for produced the thermal neutron 


fission 


Mass spectrometric ratios 


Ratio after dilution with natural 
atoms per gram plutonium 
atoms per gram plutonium 


Atoms fission per gram plutonium 


137 


9.880 
5.299 
133 
1.000 
1.000 
1.000 
1.405 
0.7193 
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TABLE 


The mass spectrometric and isotope dilution data for produced the thermal 
neutron fission 


Isotope 
142 143 

Mass spectrometric ratios 

1.000 

1.000 

Ratio after dilution with natural neodymium 

atoms per gram plutonium 2.026+0.015 1.000 

atoms per gram plutonium 1.000 
Atoms fission per gram plutonium 

1.076 


the evaluated number atoms and per gram 
plutonium each sample. From the relative yields for each the elements 
given Tables and the absolute number atoms one isotope 
each these elements given Tables the number atoms each 
fission nuclide per gram plutonium has been tabulated columns and 
Table for each the two samples. 


TABLE 


alloy 
Atomic No. atoms/g No. 
118-130 5.70* 
131 3.83 3.925 3.70 
132 1.070 5.34 5.477 5.16 
133 1.405 7.02 6.77 
134 7.59 7.775 7.32 
135 1.477 7.38 7.552 
136 1.349 6.74 6.907 6.51 
137 6.60 6.763 6.37 
138 6.43* 6.574 6.19 
139 6.05 5.92* 
140 1.110 5.55 5.61 
141 5.14 
142 0.9880 4.94 5.299 4.99 
143 0.9081 4.57 
0.7637 3.81 4.09 3.86 
145 0.6220 3.11 3.326 3.13 
146 0.5113 2.56 2.734 2.58 
147 2.056 1.94 
148 1.70 1.72 
149 0.2661 1.346 1.27 
150 1.01 1.083 1.02 
151 0.1638 0.82 0.9288 0.78 
152 0.1258 0.60 
153 
154 0.299 0.3029 0.285 
155 
156 


*Assumed yields. 


Cumulative yields based calculated number fissions 
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The yields the xenon isotopes given Fleming and Thode (1956) are 
normalized the cesium isotopes means the factor 0.925+0.010 given 
these authors for the ratio the mass region below mass 130, 
radiochemical yields have been used where possible. Other yields appearing 
Table including masses 139 and 141, which have not been obtained directly, 
have been interpolated from plot yield versus mass. 

The yield the 144 mass chain may obtained from both the cerium and 
neodymium isotope dilutions. For sample the value obtained from the 
neodymium 15% higher than from the cerium. may seen from Table 
that for this sample too much natural neodymium was used the isotope 
dilution that the ratio after isotope dilution was very nearly 
that natural neodymium. For such conditions the accuracy the deter- 
mination the number fission greatly reduced. Table 
therefore, the values for the neodymium isotopes sample have been 
obtained from the relative yields shown Table and the number atoms 
mass 144 taken from shown Table VI. For the plutonium—aluminum 
alloy sample identical yields mass 144 were obtained from both the cerium 
and neodymium isotope dilutions. This agreement obtained when value 
278-day half-life for assumed for the growth and decay corrections 
the observed mass spectrometric data both the neodymium and cerium 
given Tables and respectively. 

The yield was determined only for the alloy sample using the 
isotope dilution data shown Table The number fission atoms 
per gram plutonium given Tables and 

Two methods may used obtain absolute cumulative fission yields from 
the data given columns and Table one method the absolute 
number fissions calculated from the irradiation data and the fission cross 
section and the other method the relative yields are adjusted that 
they total 100%. The absolute number fissions the may calculated 
using the following relation (1): 

(1) 


where number atoms the nuclide, 
number fissions, 
(nv)t= integrated flux neutrons 
the flux self-shielding correction for the irradiated samples. 
equation (1), the integrated flux and the effective neutron fission 
cross section require determination. The former can related the 
disintegration rate, the cobalt monitors according equation (2): 


(2) 


number cobalt atoms, 
effective cross section for cobalt 
the flux self-shielding correction for cobalt wire. 


= 
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The effective neutron cross section, any nuclide function the 
neutron energy distribution and the variation its cross section with energy. 
The value may calculated from the following expression given 


Westcott (1958): 


where the monokinetic cross section 2200 

represents the departure the nuclide’s cross section from law terms 
its effect for Maxwellian neutron spectrum, 

epithermal index defined Westcott (1958), 

characterizes the departure the nuclide’s cross section froma law 
the epithermal region. 

For the cobalt monitors the value was taken 0.01, 36.5 barns 
(Jervis 1958), and 123° for the apparent temperature the Maxwellian 
component the neutron distribution (see calculation below). Using the values 
and taken from tables given Westcott (1958) one obtains value 
36.9 barns for the effective neutron cross section 

Using the methods previously described (by Petruska al. (1955)) 0.977 
self-shielding corrections were made for the cobalt monitors. With these the 
measured specific activities the cobalt monitors shown Table XII and the 
values the neutron flux shown Table XII have been calculated with 
equation (2). 

The value the effective fission cross section, for dependent 
the apparent neutron temperature. The apparent neutron temperature may 
calculated from the effective cross section This was obtained from 
the measured ratio the samarium isotopes the irradiated samarium 


monitor using equation (4). 


150 
where 1.86 the natural ratio the unirradiated sample, 

the observed ratio after irradiation. 

For the alloy sample, the observed 
value 123° for the neutron temperature may then estimated from 
equation (3) using 42,200 barns for (Bidinosti 1958), 0.01, and 
values and given function temperature Westcott (1958). For 
this temperature value 872 barns estimated for the effective fission cross 
section from equation (3) using 750 barns for the 2200 fission 
cross section and the values and tabulated Westcott (1958). 

Although the sample did not have samarium monitor, was 
irradiated the same position the NRX reactor and the fission cross section 
for plutonium has been assumed have the same value. With this value for 
the effective cross section the self-shielding corrections calculated 
the same manner for the cobalt monitors and the integrated fluxes given 
Table XII, the number fissions per gram plutonium shown this table 


have been calculated using equation (1). 


4 
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With the yields atoms per gram plutonium given columns and 
Table and the number fissions per gram plutonium given Table 
the yields shown columns and Table have been obtained. 

The method adjusting the absolute cumulative yields the heavy mass 
curve add 100% accomplished taking the ratios 100 times the 
number atoms per gram plutonium which are given columns and 
Table the total numbers fissions per gram plutonium which are the 

The final fission yields shown Table and represented graphically 
Fig. have been obtained this manner. 


YIELD 


FISSION 


23 127 31 135 139 143 147 1S! 155 
MASS 


Fic. mass yields thermal neutron fission 


DISCUSSION 

The accuracy the relative the isotopes given element depends 
only the accuracy the directly determined mass spectrometric measure- 
ments whereas the fission yields also depend isotope dilution measurements. 
The precision the relative vields given Tables are each case 
better than and agree within these limits with the values previously 
published (Wiles and Tomlinson 1956; Krizhanskii 1957) except three 
mass regions. The most serious difference the ratio shown 
Table Since the stable isotope which occurs nature may con- 
sidered that high values this ratio are the result stable contamination. 
the previous determination from this laboratory (columr Table I), only 
partial recovery the fission products was achieved. Recent unpublished work 
this laboratory has established that the mobility the 5-day precursor 
can materially alter fission ratios cesium isotopes oxide 


samples. 
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TABLE 


Cumulative yields and other literature values 


yield 


Atomic Steinberg al. 
number alloy Average Wiles (1956) (1955) 
118-130 5.70* 5.38 

131 3.80 3.73 3.6 
132 5.30 5.21 5.26 3.79 4.9 
133 6.96 6.83 6.90 4.97 5.0 
134 7.52 7.39 7.46 5.37 

135 7.32 5.22 5.5 
136 6.69 6.56 6.62 4.77 

137 6.54 6.42 6.48 4.94 5.8 
138 6.25 6.31 5.38 

139 5.97* 6.61 5.4 
140 5.50 5.66 5.58 7.36 5.36 
141 5.10* 6.94 4.9 
142 4.90 4.97 6.62 

4.50 4.61 5.98 
144 3.78 3.89 3.84 5.00 
145 3.08 3.16 3.12 4.07 

146 2.53 2.60 2.57 3.36 

147 2.02 1.96 1.99 2.81 

148 1.69 1.73 2.27 

149 1.32 1.28 1.30 1.81 

150 1.00 1.03 1.02 1.31 

151 0.814 0.790 

152 0.625 0.606 0.616 0.88 

154 0.297 0.289 0.293 0.40 

Total 100% 100% 100% 100% 


*Assumed yields. 
Note: Summation of measured yields together with values interpolated from smooth mass-yield curve by 
Steinberg et al. (1955) totals 93%. 


TABLE 


Number fissions the plutonium samples 


Specific activity Number 

cobalt monitors, Self-shielding 


*Average for two monitors. 


The value obtained Krizhanskii a/. (1957) can not assessed without 
further information, but the new value presented Table has been con- 
firmed with measurements from several entirely different irradiations. 
similar type variation, but lesser extent, has also been observed 
mass 140 (see Table which apparently does not depend 


contamination. 
further difference observed the samarium yields masses 152 and 
154 (see Table The vields given Krizhanskii (1957) were corrected 
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assuming contamination arising from natural samarium. fact the 
with and hence this correction should not have been made. The values 
Krizhanskii al. (1957) without correction for correspond closely the 
values obtained this work. 

Fritze and co-workers (1958) their simultaneous determinations xenon 
and krypton from the same alloy sample obtain cumulative fission 
yield 3.79+0.11 for This good agreement with the value 3.77 
obtained this work from the measured yield and the relative yields 
and previously given Fleming and Thode (1956). 

The final values the cumulative fission yields given Table 
have been normalized that the yields total 100%. The values given 
Table depend the calculation the total number fissions occurring 
each sample. the case the sample, where the neutron temperature 
was not measured (hence the fission cross section doubtful) and where the 
self-shielding correction amounted about considered fortuitous 
that the vields shown Table total For this sample the values 
shown Table are therefore preferred. the case the alloy 
sample these factors were more carefully controlled and therefore the 99.11% 
the vields Table could well measure the cumulative 
errors the several determinations the values the few yields which have 
been assumed. The final values for this sample shown Table have never- 
theless been normalized 100% since they are then independent the 
calculated number fissions, factor which dependent the relative cross 
sections and and the self-shielding correction. The agreement 
with the yields shown Table taken measure the absence loss 
fractionation fission products from the sample. 

The yields Wiles a/. (1956), also shown Table XI, are now assumed 
error from fractionation which occurred the partial extraction the 
fission products from the The radiochemical yields summarized 
Steinberg a/. (1955) are difficult compare with the present values since 
their summation only 93%. 

The values the yields obtained from the two independent determinations 
shown columns and Table agree better than most masses 
and the average values shown column Table are considered have 
absolute accuracy better than each mass. 
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THE SUPERCONDUCTING TRANSITION 
POLYCRYSTALLINE 


ABSTRACT 


The influence grain size the breadth the superconducting transition 
polycrystalline tin has been studied. The transition has been detected using 
a-c. induction method. range grain sizes from 0.2 was used and 
the corresponding range transition breadth from millidegrees. The 
influence the grain boundaries broaden all the transition curves towards 
higher temperatures from well-defined end point which common all the 
curves. 


INTRODUCTION 


view the interest attached the influence lattice imperfections 
the superconducting transition and the apparent lack information the role 
grain boundaries specifically, experiments have been done study the effect 
grain size the superconducting transition. Work thin films has shown 
marked dependence the superconducting properties the lattice structure 
but has not apparently isolated the influence grain boundaries alone. The 
only work done hitherto bulk specimens seems that Haas and 
Voogd (1931), who quote results for only three samples and, these, only two 
are related purity. figures for the actual grain size are given. The problem 
the influence grain boundaries the superconducting transition 
interest, not only for its own sake considering the mechanism the transi- 
tion but also aiding the evaluation experiments other properties 
superconductors which are performed polycrystalline materials. The aim 
the present work, therefore, repeat and extend the measurements 
Haas and Voogd attempt clarify the nature the transition 
broadening found polycrystals. 


SPECIMENS 

The specimens were made from zone-refined tin obtained from the Depart- 
ment Physical Metallurgy, University Toronto, and two samples taken 
from the same bar were used. One was cut from the original bar 
composed almost entirely single crystal with only one longitudinal grain 
boundary. The other was cut random and severely worked hammering 
achieve the smallest possible grain size. Each was turned into cylinder 
approximately diameter and long. The worked surface produced 
machining was etched away, thereby simultaneously revealing the grain 
structure for grain-size measurements. Each specimen was annealed 100° 
for hours ensure that the lattice within the grains was free stress 
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possible. After measurements its superconducting transition described 
below this second specimen was annealed temperature just below the 
melting temperature for the time required increase the average grain 
diameter. The transition experiments procedure were 
repeated the specimen until further grain growth could achieved. 


METHOD 

The transitions were detected magnetically method described 
Webber, Reynolds, and McGuire (1949). The specimen was mounted vertically 
liquid helium bath and was enclosed close-fitting secondary coil. 
oscillating longitudinal magnetic field was produced solenoid wound round 
the outside the helium Dewar. This primary coil gave about gauss per 
ampere and had additional compensating turns the ends ensure constancy 
the magnetic field over the specimen volume within .2%. This coil was 
driven audio oscillator and the value the magnetic field was calculated 
from current readings. The signal from the secondary coil was led the top 
the cryostat coaxial cable and fed amplifier which could fitted with 
narrow band-pass filters. The output the amplifier was read RMS 
vacuum tube voltmeter. 

The electromotive force (e.m.f.) induced the secondary coil depends the 
state the specimen which forms the core. the normal state the flux linked 
that due the oscillating field diminished only the diamagnetic influence 
the eddy currents induced the specimen. the purely superconducting 
case the flux linked should very close zero but small contribution remains 
due leakage flux through the small space between the coil and the specimen 
and through the penetration depth the magnetic field into the specimen. The 
main difficulty with this method that, with specimen high purity, the 
flux linked the normal case may diminished eddy current shielding 
that the change over the superconducting transition may quite small. 

The conditions the intermediate state are much more complicated. Since 
these intermediate state effects influence the transition which observed 
this technique, and since the method being used not merely detect super- 
conductivity but study the transition itself, interesting analyze this 
situation order know precisely what being measured during the transi- 
tion. The process is, the general case, one repeated magnetization the 
superconducting state, passage through the intermediate state, and destruction 
superconductivity alternate directions. This illustrated Fig. which 
the peak value the external field shown relation portion the 
diagram the specimen. the external field varies time during one 
cycle from the point the specimen remains superconducting until certain 
fraction the value the critical field that temperature reached 
From until the full value reached the specimen the intermediate 
state after which fully normal until the value reached again with 
diminishing external field. During this whole process there has been certain 
value giving the signal which observed. The observed transition 
then consists the variation the fundamental component the 
temperature lowered from until the peak external field 
longer able generate the intermediate state. 


— 
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He 


NORMAL 


INTERMEDIATE 
STATE 


Cc A Te T 


Fic. Comparison the applied oscillating magnetic field with the critical field curve. 


The values and can calculated for any portion this transition 
region the magnetization during the superconducting and intermediate states 
known. Shoenberg (1937) has measured the magnetization cylinder 
ratio unity longitudinal magnetic field through the super- 
conducting and intermediate states. the case non-ideally shaped body 
the magnetization phenomena are complex and largely governed the 
behavior trapped flux. These effects depend considerably the shape and 
condition the specimen but still interest use magnetiza- 
tion curves comparison with the present results even the cylinders are 
different ratios. 

Shoenberg’s curves intensity magnetization can replotted give 
function and differentiation gives the required variation 
with external field. The form these curves depends the relative values 
the peak external field and the field required initiate the intermediate 
state. Two them, one complete cycle, for and for are 
given Fig. Fig. are pictures taken oscilloscope the amplified 


-2- 
He = He 
TIME 


Fic. Predicted variation induced e.m.f. the secondary coil two temperatures 
during the transition. 
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Fic. the secondary coil output during the transition. They are numbered 
order decreasing temperature. Frequency 1000 cycles/second. gauss. 


output from the secondary coil (without filters). They were taken intervals 
through the transition from point just after the transition 
until just before the signal becomes its residual constant value and represent 
the output for various values the ratio These observations were made 
the largest grain-size specimen using comparatively large oscillating 
about gauss, illustrate the intermediate state properties. The 
predicted curves Fig. should compared with pictures and Fig. 
which refer approximately the ratios quoted. Two points interest arise. 
The first that, anticipated, the trapped flux region which dominates the 
=H, curve only approximately represented the predicted curve. One 
the possible contributions this discrepancy suggested the second 
feature the curves which exhibited the curve for The predicted 


as 
6 
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curve for this ratio has distorted region arising from the passage through the 
superconducting and intermediate state regions when the external field 
passing through the region near zero, i.e., maximum. However, the 
observed curve shows kink displaced from the maximum the e.m.f. 
presumably ascribed lagging the flux changes arising from finite phase 
boundary velocities. This supported comparison the present observa- 
tions frequency 1000 cycles per second and the similar oscillograms 
shown Galkin and Bezuglyi (1955) for frequency 200 cycles per second 
where the displacement less. The signal which was measured during the 
transition the fundamental component the curves but attempt 
has been made analyze the curves for comparison with the measured signals. 
This the nature the transition which observed this method and 
used study the superconducting transition itself and measure 
transition widths there the obvious necessity work with the smallest 
possible values the present work values down 0.1 gauss could 
used and this would correspond intermediate state width only 
millidegrees, which small compared with the observed transition widths 
about that this instrumental transition 
width not influencing the observed transition widths comes from observation 
the transition with various values H,. found that the observed transi- 
tion width can reduced reducing the value only certain point. 
Thereafter the transition width independent the primary field and all 
results illustrated refer this limiting condition. any case, attempt 
provide shielding from the earth’s magnetic field was made and this could 
account for transition width about millidegrees. Even this, however, 
not serious factor the case the small crystal specimen although 
presumably does contribute more significantly the transition for the largest 
grain-size specimen. 


PROCEDURE 


The transitions were studied measuring the output the secondary coil 
while slowly reducing the temperature. The pressure was read mercury 
manometer which measured absolute pressures and which could read 
1/10 mm. Corrections the temperature for hydrostatic pressure head 
helium above the specimen were made and temperatures were calculated 
using the 755, temperature scale (Clement al. 1955). The transitions were 
repeated over range values and for the two frequencies 400 


and 1000 cycles per second. 


RESULTS 


The transition curves for the various grain sizes are shown Fig. The 
quantity plotted vertically the secondary coil output normalized that the 
drop signal over the transition from unity zero. Using this normalization 
dependence the transition frequency could detected and all the 
curves refer values the peak external field gauss) low enough 
negligible affecting the transition width. The curves are composed 
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FLUX LINKAGE 


3.70 3.72 3.74 3.76 378 
TEMPERATURE 


Fic. for the various grain sizes. The average grain sizes are, reading from 
left right: approximately single crystals; 5.5 mm; (c) 3.0 mm; (d) 1.4 mm; (e) 0.2 mm. 


numbers points ranging from for the sharpest transition about 100 for 
the widest obtained from repeated runs. The scatter the points never more 
than total range millidegrees temperature and .04 susceptibility. 
These extreme values are exceptional and the normal scatter about half the 
value quoted. The influence the transition the grain-size width illus- 


TRANSITION BREADTH 


2 4 6 8 10 (2 MM 
GRAIN SIZE 


Fic. Dependence total transition width grain size. 


trated Fig. this the transition width the total width the transition 
from the point which the onset could just detected the point which 
the induced e.m.f. had just become constant (both detectable about one part 
one thousand the total change). The uncertainties transition width are 


10 
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estimated from the scatter the measured points the curves and the un- 
certainties grain size from the observed range diameters. The specimen 
which was almost single crystal difficult define this way and the value 
quoted for the mean grain size estimate the range single crystal 
current path for the circumferential currents induced the experimental 
procedure. 

Another specimen which was much more nearly single crystal (containing 
only small amount striation substructure) gave total transition breadth 
only millidegrees. However, came from different batch material and 
not quoted part the series illustrating the influence grain boundaries 


alone. 


DISCUSSION 


The general result the experiment confirm the original observation 
Haas and Voogd that the influence grain boundaries broaden the 
transition only the direction increasing temperature away from fairly 
well-defined vanishing point resistance. Although not visible the 
scale Fig. all the the superconducting state within about 
millidegrees between 3.697 and 3.699° The broadest transition reported 
Haas and Voogd for wire covers about millidegrees, while the broadest 
one presently reported for bulk specimen about millidegrees. This last 
probably close limiting value the specimen had about the smallest grain 
size which can obtained pure tin. 

For single crystal Webber, McGuire, and Reynolds report transition 
breadth millidegrees and this can compared with the value 
millidegrees for the largest grain size the series. stated above, specimen 
freed from the last traces large angle boundaries and most the striation 
substructure gave total width about millidegrees that the first intro- 
duction misorientation has relatively large effect the transition breadth. 

The problem thus explain the nature the partial superconductivity 
observed above the bulk transition temperature. not yet possible say 
whether this effect taking place the grain boundaries themselves 
whether can arise from thermal contraction mismatch between areas 
crystal misorientation. For ideal case misorientation Faber (1957) has 
evaluated the transition temperature shift tin the order 100 
millidegrees upwards. This certainly covers the observed magnitude the 
effect but the extension the treatment the genuinely polycrystalline case 


not easy. 
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THE SPHERICAL OSCILLATOR 


ABSTRACT 


The degree consistency oscillator model nucleus examined 
means type calculation based simple form inter- 
nucleon potential valid low energies. effective mass equal 0.757 times the 
mass free nucleon used. The oscillator wave functions are found not 
far from self-consistent and the oscillator frequency derived physically reason- 
able, but the bound the binding energy not good. also shown that the 
oscillator wave functions are good approximation for the state functions 
particles bound finite potential well having the shape cutoff oscillator 
that the Hartree-Fock calculation can used prescribe shell model 
potential. 


SECTION THE PROBLEM 

Frequently problems nuclear physics necessary assume some 
radial wave function represent the motion single nucleon nucleus. 
For such cases oscillator function often chosen (e.g. Talmi 1952; Tauber 
and 1957) principally because its properties are well known and its 
analytic behavior relatively simple. these cases one concerned obtain 
estimate some matrix element and not committed any statement 
about the potential which the nucleons move. 

However, there have been several more detailed studies which properties 
physical nuclei have been related oscillator shell model. The earliest 
shell model level orderings obtained with coupling (Mayer 1949; Haxel, 
Jensen, and Suess 1949) used that model. More recently Nilsson (1955) has 
made extensive calculations levels spheroidal oscillator wells and these 
have been used combination with the collective model Bohr and Mottelson 
(1953) for the prediction spins and parities nuclear ground states (Mottel- 
son and Nilsson 1955) and for quite detailed comparison with excited states 
light nuclei (for example Litherland, Paul, Bartholomew, and Gove 1956), 
but there not, the literature, any real justification for the use oscillator 
wave functions for nucleons nuclei. 

shall try show that there is, fact, some justification for using the 
oscillator form the central potential single-particle model nucleus. 
The method determine, quantitative fashion seems easily 
possible, the departure the oscillator potential from self-consistency 
calculation using simple non-singular internucleon potential 
(equation 17) which correctly describes the two-nucleon interaction low 
energies, i.e. yields the description the deuteron and various its inter- 
actions and Schwinger 1951) and fits the low-energy proton—proton 
differential scattering cross section (Hall and Powell 1953). 
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Our purpose modest; such calculation can not expected correctly 
determine, for example, the binding energy nuclei the true density 
nuclear matter. This requires much more refined methods, including the 
correlation nucleon motions and the proper description high-energy inter- 
actions within the nucleus. The procedure not applicable 
the singular potentials computed Gammel and Thaler (1957) and would 
very difficult apply the complicated potentials Signell and 
Marshak (1958). 

the other hand, for examination the most suitable wave functions 
use the unperturbed shell model one may permitted use simple 
methods and potential fitted only low energies. Previous calculations 
this general nature have been made Rotenberg (1955) and Talman 
(1956). Both these calculations started with the wave functions square 
well and then improved the function and potential interaction. both 
cases, but particularly Talman’s calculation, the resultant potential appears, 
superficially least, more like cutoff oscillator than like the original square 
well. 

The calculation presented Section was designed examine the useful- 
ness the oscillator functions initial approximation and see what 
results can obtained that approximation the use momentum- 
dependent shell model potential, i.e. effective mass less than unity for 
nucleon nucleus (Brueckner and Gammel 1958). The calculation does not 
contain iterative process. Instead have computed matrix elements be- 
tween single-particle-occupied states the Hamiltonian opera- 
tor, taking exchange completely into account. These elements can then 
compared directly with the corresponding elements the single-particle 
Hamiltonian. The comparison shows qualitatively that the oscillator wave 
functions are nearly self-consistent. absolute result obtained, because 
the nature the calculation does not justify the time and expense attempt- 
ing the self-consistency calculation. 

However, one does obtain explicit statement, known unperturbed 
oscillator shell model, the energy level scheme appropriate inter- 
nucleon potential valid nuclear energies. This was not provided the two 
earlier computations and may interesting compare with the calcula- 
tions binding finite nucleus. After the work reported here was com- 
pleted Brueckner (1959) announced self-consistent calculation using 
the effective two-body potential derived from his treatment infinite nuclear 
matter. The comparison his results with ours gives indication the 
effect properly treating the high-energy interactions. 

clear that any reasonable shell model potential must zero outside 
the nucleus. Therefore first establish (in Section that the bound eigen- 
values and eigenfunctions potential are negligibly 
different from those the corresponding states complete oscillator. 
Once this result established will use simple oscillator functions since 
only bound states are treated this paper. 
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SECTION THE CUTOFF OSCILLATOR POTENTIAL 


preliminary, will here examine the single-particle bound state wave 
functions which satisfy the Schrédinger equation 


V=0 


(2) represents the true mass nucleon and the effective mass 
defined according the present custom nuclear physics, i.e. that the 
approximation the true total Hamiltonian given shall con- 
tain momentum-dependent term 


The mass parameter for the single-particle wave equations then becomes 

defined (2d). and are coupling constants used Nilsson (1955). 
put 


function spin and angle variables with the properties indicated 
the subscripts; 


(4) X ijm = (I, 3 2) m—v, vijm) (Ob) x4 


(4) the bracket representation used for the coefficients; 
normalized spherical function and spin function. The differential 
equation for 


where 


The solution (5a), which well behaved the confluent hyper- 
geometric function 


1 ( 1 
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the notation Bateman (1953, Vol. 248). Thus the bound states 
equation (1) are 


(7b) = Cu Kisy (xr) X jm (9, in r R,. 


(76) that Bessel function imaginary argument defined 
Watson (1952) which decays exponentially for large and 


(7c) 


and are normalizing constants. 

The energy eigenvalues the bound states are determined the require- 
ment that the logarithmic derivative continuous R,. This 
condition, terms ak?, may written 


where all functions are evaluated R,) and the prime represents 
differentiation with respect the argument the function which 
attached. 

define 


(9) 
that 
and define also 


then the term the left (8) may evaluated function from 
the recursive relation 


So(v, 
Si(v, 


Finally there the relation (Bateman 1953, Vol. 254) 


(a, x) = c ad(atl, c+1;x) 


that (8) becomes equation for 


For given and each solution (12) defines the radial 
the potential (2). increased the values 
approach i.e. the energy eigenvalues isotropic oscillator; 


= 
= 
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terms the parameters tends toward negative integer, 
Hence the effect the cutoff may represented the value n+a 
and also the energy shift. 

Table displays the parameters which define bound states cutoff oscil- 
lator well when 0.6089 and for the indicated values and 


TABLE 
Bound energy levels the cutoff oscillator potential equation (2) with 0.6089; 
—n+e. the increase the energy the cutoff oscillator over 
that the complete oscillator 
11/2 


| 
| 
| 
| 
| 
| 
| 
| 
| 


The parameters Vo, and have been chosen that comparisons can 
made with the Hartree-Fock values Section For simplicity has 
been given the same value for all shells. This case completely typical; 
unnecessary give tables for other values Vo/hw. 

may seen that there remarkably small change energy due 
cutting off the oscillator potential, and this result not modified the 
presence terms Even levels lying very near the top the well 
are but slightly changed. The potential (2) would appropriate for the 
neutrons model nucleus. The proton potential should include Coulomb 
barrier and the shifts proton levels would even smaller than those 
shown Table 

interesting observe that the direction the energy shift varies 
with the magnitude easily seen that when all bound 
levels are lowered energy, but decreases one finds that some levels 
are raised. This caused, course, the changed condition the slope 
the radial function R,. the logarithmic derivative imposed the 
Hankel function the external region smaller magnitude than that 
the original oscillator function, the energy the state decreased, and vice 
versa. Since the effective mass ratio unity outside the nucleus the imposed 
slope independent while the slope the oscillator function depends 
the effective mass used the internal region. 

Not only are the bound energy levels very nearly undisturbed the cutoff, 
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but also the wave functions are very little changed. This may demonstrated 
comparing the positions the nodes the wave functions the cutoff 
and complete oscillator potentials, i.e. have consider the roots 

this equation has just real positive roots when (Bateman 1953, 


288). 
This study can aided consideration the polynomials 
and 


n k 


k=0 


tween the values the corresponding terms and Conse- 
quently the values the polynomials are similarly intermediate, and also 
are the positions the roots. Thus denote the root 


Now the difference 
denote the root (13), then 
(146) Vn—2x €) < Xn—2%(, €) 
€) > €). 

The equality sign holds when that from and 
have 
(15d) 0) < 


When becomes unity becomes infinite, i.e. coalesces with the 
singularity infinity, while 


€). 


From (15) see that 
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Thus, although this discussion has been limited positive values the 
analyticity requires that (16) hold also for near zero 
roots, while there are just roots 0). this case represents 
infinity and relation (16) does not follow from (15) This root, 
however, lies the region and not node our wave functions 

For the values shown Table see that the radial zeros are 
shifted about 0.5% (about change from the zeros the com- 
plete oscillator functions. The principal difference between the functions (7) 
and oscillator functions occurs the regionr (7) describes 
surface, but the distribution particles within the nucleus very little 
different from the distribution defined oscillator functions. 


SECTION HARTREE-FOCK CALCULATION 


this section describe (denoted H.F. sequel) calcu- 
lation the ground state model nucleus. assume nuclear wave func- 
tion constructed from one-particle oscillator functions. the last section 
have learned that these are negligibly different from the corresponding single- 
particle states cutoff oscillator potential, that they form reasonable 
basis for the calculation. Moreover, well known that nuclear model 
based these functions displays many the properties physical nucleus, 
particularly when combined with the collective motions the theory Bohr 
and Mottelson (1953) when treated the methods Elliott (1958). 

will compute the matrix elements between occupied levels the H.F. 
Hamiltonian, derived from the two-body potential, 


(17a) V(r) 


(17) the Majorana operator which exchanges position co-ordinates. 
The constants are those used Lane, Thomas, and Wigner (1955). 

The potential (17) not suitable nuclear potential all energies. Recent 
high-energy and scattering experiments have 
best been described internucleon potential containing repulsive core 
and much more complicated form outside the core. non-integrable core 
(Gammel and Thaler 1957) normal type calculation 
leads divergences, and the potential Signell and Marshak (1958) too 


complicated for easy analysis. 

One can give qualitative argument for using H.F. calculation the 
simple form (17) together with momentum dependence; the basis the 
argument that the true potential first replaced pseudo potential, and 
the latter used for the H.F. process. Pseudo potentials have been 
point operators (Huang and Yang 1957; Lee, Huang, and Yang 1957) for 


| 
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hard spheres, but different kind pseudo potential can imagined. For 
unbound spinless particles, the specification all phase shifts function 
match the effect given potential all unbound particle motions, 
suffices match all asymptotic phase shifts. If, addition, the given potential 
vanishes on, can exactly treated outside, some surface suffices 
match the phase shifts that surface. This may done either specifying 
boundary condition there, specifying some other momentum-dependent 
pseudo potential within the boundary which will produce the same set 
phase shifts. example such pseudo potential for very simple original 
potential given appendix this paper. 

This computation would become more complicated were attempted for 
the recent internucleon potentials. The presence tensor and spin-orbit 
forces requires the matching scattering matrices instead simple phase 
shifts, and the choice suitable boundary surface not easy. not 
here intend attempt such calculation; the concept introduced merely 
point out that computation with the simple form (17) not wholly with- 
out sense. That form, combined with tensor term, was proposed Fesh- 
bach and Schwinger (1951), and fitted the low-energy data the two- 
nucleon system. Hall and Powell (1953) fitted their scattering 
data Mev using the same form and determining the remaining free 
parameters. Thus the potential, which (17) the 
central part, may regarded experimentally determined pseudo poten- 
tial valid low energies. The tensor part the potential 
would make contribution the energy system with saturated spins, 
that the filled shell systems which discuss its effect would always 
small, i.e. would have (small) effect only the states incomplete 
Hence legitimate our calculation use only the central term, 
i.e. equation (17). 

The concept pseudo potential outlined above not very different 
principle from the more elegant and more rigorous concepts 
theory (see for example Bethe 1956; Gomes, Walecka, and Weisskopf 1958). 
That theory has been used show that infinite nuclear matter complex 
internucleon potential approximately equivalent one-body potential 
which includes term proportional (Brueckner and Gammel 1958). 

shall not here attempt justify the use any particular value the 
momentum dependence the one-body potential. The fact that effective 
mass less than unity should used shell model nucleus now 
well established; method which, principle, should permit calcu- 
lated within the framework the Hartree-Fock process has been proposed. 
treat the ratio (cf. equation parameter the computation. 

remains briefly describe the computation the H.F. potential matrix. 
shall calculate only for systems which both neutrons and protons form 
closed coupled angular momentum shells. The trial wave function for the 
system particles then Slater determinant composed from single-particle 


wave functions, 
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where 


The function (18) differs from (38) that here use the wave functions 
the complete, not the cutoff, oscillator and (18) also includes function 
the isotopic spin order that the total nuclear wave function may 
written single Slater determinant. Since the potential (17) charge 
Laguerre polynomial the notation Bateman (1953, Vol. 188). The 
function angle and spin has been defined Section The coupling 
imposed, say, the Nilsson (1955) coefficients, but this not significant 
for our H.F. calculation. The justification the spin-orbit coupling term 
shell theory unsolved problem long standing which the present 
calculations make contribution; since limit ourselves closed shell 
systems only the radial dependence the central potential can discussed. 

Now let stand for the five co-ordinates one particle. set up, 
and shall seek minimize, the integral 


which represents the Slater determinant formed for neutrons and 
protons (V+Z .1) from the single-particle functions (18). the 
kinetic energy operator and the two-body potential that defined 
(17). The integral includes spatial integration and sums over all 
spin and isotopic spin variables. The reduction (19) when spin and orbital 
angular momentum are uncoupled obvious; essentially the same reduction 
possible for our coupled wave functions. The results are particularly 
simple when includes the Serber exchange mixture contained 
(17). note first that the potential terms (19), after expansion the 
(normalized) Slater determinants, use the symmetry the potential sum, 
summing over isotopic spin, and application the exchange operator 
become 


(20) represents the three spatial co-ordinates particle the states 
longer contain the factor and the sums over extend 
over all occupied states. now consider that part the sums (20), ex- 
tending over one two closed shells, which the radial part the wave 
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function constant. The angle and spin function written out 
equation Section (2). have 


(21 ) = X Ce) 


=m 


obtain (21) have used the normalization property and the symmetry 
under interchange (jm) with the coefficients. 
Equation (21) provides the means for reduction the first two sums 
(20). The other two sums arise only when the states are filled nucleons 
like charge. the ground state will generally true that least one 
nucleons the same kind. that case the double sum can simplified. 
Suppose and are both occupied (this procedure also 
applicable when only one shell occurs); then 

Equation (22) follows from the ortho-normal properties the 
coefficients. Since may write, 
for the case when least one completely filled, analogy with (21) 


mm’ 


aq 


Note that (23) does not contain sum over 
shall use (21) and (23) together with the addition law 


q 
simplify (20). the system has two more unpaired shells 
filled, unfilled, containing nucleons like charge, then (23) 
approximation for those particular shells. 
Putting 


(25) Zn iq(X) = Riilr) Yi(4, >), 


have 
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Except for the occupation numbers equation (26) 
now the same one would have obtained from Slater determinant func- 
tions with uncoupled spin and orbital angular momentum. This the only 
remaining effect the coupling and the matrix elements 
between states Z,, are course independent the coupling 
constants used. The subscript represents the triple quantum numbers 
The sum contains repeated terms for multiply occupied states. 

The H.F.-linked system equations obtained requiring that 
the integral (26) stationary under variation the single-particle wave 
functions This standard procedure which leads system 
equations 


Equation (27) can, principle, solved for the oscillator 
nucleus were self-consistent the H.F. sense, the which satisfy (27) 
would turn out just oscillator wave functions, and the spherically 
symmetrical systems under study this solution independent the 
and orbital momentum dependence the potential used define the original 
functions. 

The direct terms contain factor 6,, and can easily compared 
function with the original oscillator potential. They indeed show 
approximately quadratic dependence out the nuclear boundary with 
smoothed cutoff zero potential energy there. The direct potential is, how- 
ever, much shallower than the original oscillator form. 

order obtain fairly simply definite comparison between the com- 
plete H.F. equations (27) and the original single-particle equations, have 
decided compute the (single particle) matrix elements (27) and 
see how they differ from the matrix elements the original isotropic 
oscillator. This simpler calculation, and more definite than the usual 
comparison the average over V,,(x) with the oscillator potential. 

The matrix the kinetic energy operator well known, have only 
compute the potential matrix whose elements are 


g 


Now the spatial dependence the potential (17) may written (Mott and 
Sneddon 1948, 387) 


A 


which and these are some order; the angle between 


and X», and are Bessel functions imaginary argument 
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defined Watson (1952). When (29) substituted (28) obtain the 
following sums integrals over angle: 


qo 
= P,(2\+1) 6001, 1,59,¢,(2le +1); 


qo 


The matrix element (28) then reduces sum integrals over the two 
radii, and continue the separation into direct and exchange terms. 


The integrals (31) are sums Slater integrals the form 


0 
those obtained replacing sinh cosh ur. 

The calculation was performed the Chalk River datatron 
evaluating integrals form for sufficient range values and 
and then finding the combinations required for the sums The Coulomb 
interaction between protons contains Majorana operator obtained 
from the same formulae setting replacing —e? the direct 
term, and +e? the exchange term; and limiting the sum proton-occupied 
states. The integrals (82) were evaluated integration. 
This relatively simple since they are essentially single integrals. 
results were checked comparison with analytic values special cases and 
recursion relations. For small values and the computed 
integrals were accurate six significant figures; for larger values this 
improved about one (Simple programs the datatron are limited 


eight figures.) Since these integrals have combined with varying 
coefficients and signs the matrix elements are probably accurate only 


Us 
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three four figures, and fact the exchange potential matrix elements in- 
machine values these elements fluctuate widely with small changes 
showing that significance the result almost completely lost. 


SECTION THE CHOICE PARAMETERS 

order carry out the calculation the matrices (31) necessary 
specify the oscillator parameter indeed one may choose different 
values and for the neutron and proton wave functions; they are para- 
meters the initial one-particle shell model wave functions, now con- 
sider that model. 

The cutoff oscillator potential Section defined three parameters, 
the central depth Vo, the cutoff radius and the effective mass parameter 
any reasonable nuclear model based this kind potential should 
expect that the central depth the neutron potential is, first approxima- 
tion, the same all nuclei, and that the cutoff radius given form 


(33) 


not clear that the momentum dependence the potential defined 
the ratio should also the same all nuclei, but implied 
the two previous assumptions and the systematics nuclei. From the con- 
dition that the central potential for neutrons vanish obtain 


(34) 


The average binding neutron, given the sum the expectation values 
the kinetic plus half the potential energies the occupied states is, 
system containing neutrons, 


which the number the oscillator level, and the and 


the notation Nilsson (1955). The sum runs over all neutron-occupied 
N 


states. have shown that can, with very little error, choose the 
the integers which define oscillator levels. that case when oscillator 
2n+/ among the occupied states, 


and 

N 


Then have, neglecting the spin and orbital terms, from (33), (34), 
(35), (36) 
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Now very nearly constant all nuclei. Since see from (36a) that 
(K+2) varies like the required constancy obtained only 
least very nearly the same all nuclei. choose define 
the same all nuclei that the model potential for neutrons any 
nucleus can specified three parameters a,, Vo, and 

The model potential for protons can easily found making suitable 
allowance for the Coulomb forces. the charge density supposed de- 
crease the square the radius within R,, then the Coulomb potential 
This cavalier treatment can continued replacing the quartic 
the quadratic form 


which can added the neutron potential form proton potential with 
central depth 


and with level spacing 


When have chosen set parameters for one nucleus the model described 
above allows extend the choice all nuclei. 

The nucleus Zr® was chosen for initial study. large enough for the 
model concept valid but not large that the matrix elements (31) lose 
too much accuracy. After some preliminary investigation the matrix the 
H.F. Hamiltonian for Zr® was evaluated 0.20, (.02) 0.26. 

This range values was indicated nuclear dimensions. From the 
known relation 


where the principal quantum number the spherical oscillator 
function can determine the r.m.s. radius the nucleons the oscil- 
lator nucleus. For the neutrons this ranges from 


0.886; 1.144 ata .26. 


Here the radius the sphere uniform density having the given 
r.m.s. value radius. The parameter used Fregeau (1956) fit Stanford 


data electron scattering from corresponds for .210 although 
later fitting (Ehrenberg 1959) corresponds .182 for Zr. 


general the H.F. matrix elements varied linearly with such wise 
that the mean binding per neutron varied from —2.90 Mev 0.20) 
—5.64 Mev .26). The element departs from this behavior 
some extent; decreases more rapidly for .22 than the rest the 
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elements, i.e. the large values the almost uniform spacing diagonal 
elements (cf. Tables and slightly disturbed the highest level, and 
the total binding energy increases somewhat more rapidly with increasing 
H.F. matrix ‘‘most nearly model nucleus can 
specified. The interval between the and diagonal elements the H.F. 
Hamiltonian defined 2hw. then determined from the diagonal 
element. The other model parameters are then also determined, for 


and 


1, 
(2) /h?)?/aA 3, 


The matrix the model Hamiltonian can then compared with the 
matrix the H.F. Hamiltonian. This comparison also yielded definite 
indication the best value use, but because the above-mentioned 
slightly non-uniform behavior the element the spacing diagonal 
elements was more nearly uniform .224 than other values the 
range studied. 

The range values approximately the criterion nuclear 
density; within that range the mean binding energy increases with increasing 
(decreasing nuclear size) but, otherwise, the results are essentially indepen- 
dent would have been more satisfactory determine minimizing 
the energy some criterion consistency. However, such criteria are not 
over the reasonable limits variation nuclear density; so, lacking 
any better indication have decided use 0.224 basis for sample 
computation. The value the effective mass alters very little (5% from 
.26) and the “near the model does not change 
significantly. Since these are the points interest the calculation the par- 
ticular value chosen not important. 

The model nucleus from which values and are drawn defined 


1.603 fermis, 
88.22 Mev, 
u*/M 0.6089 (M*/M 0.7569). 


SECTION DISCUSSION 

The values calculated for the parameters (39) are presented 
Tables and for the four nuclear examples choose 
again example for more detailed discussion. 

Table shows the machine output for this case using the potential 
(17) including the constants defined there. The complete H.F. Hamiltonian 
for also shown and this matrix for the other nuclei given Table 
These matrices contain only the elements between single-particle states 


which are occupied when the model nucleus unexcited. The complete solu- 
tion the H.F. equations would require the matrix elements between all 
members the (complete) set oscillator functions. Because the spherical 


The non-zero matrix elements for Zr® Mev. The rows and columns are defined 


BLE 
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the quantum numbers the radial functions, which are defined the model para- 


meters 1.603, 88.22 Mev, u*/M .6089 


Neutrons 
01) 
(04| 04) 
Protons 
(03 


47.66 —72.30 6.97 —65.33 
—20.52 —40.54 —61.06 11.62 
—16.48 —33.98 —50.46 16.26 —34.20 
—16.24 —34.83 —51.07 16.26 —34.81 
—12.62 28.03 —40.65 20.91 
—12.04 —41.60 20.91 —20.69 
8.36 —22.76 —31.12 25.56 5.56 
1.680 8.721 —10.401 5.691 4.710 
2.405 —12.773 7.347 5.426 
Model 15.26 
23.12 —50.99 11.15 
—18.68 —22.11 —40.79 15.60 —25.19 
—18.19 22.64 —40.83 15.60 25.23 
—14.00 —16.79 —30.79 20.06 —10.73 
—17.90 —30.89 20.06 —10.83 
— 1.723 — 7.490 — 9.213 5.460 — 3.753 
2.621 9.112 —11.733 7.049 4.684 
Model hwp 14.64 


TABLE III 


The non-zero elements the H.F. Hamiltonian matrices, Mev, for some other 
The numerical values the exchange contributions the elements marked (?) are doubtful 
due truncation errors the machine calculation. The model constants are 1.603, 


Model 


—54.65 —44.25 —36.19 
—40.94 —31.14 —19.43 
—28.71 18.46 
2.42(?) 
6.22 
—16.65 
6.66(?) 
0.280 
4.669 
5.328 
6.106 


12.58 


—49.19 
—10.75 
5.152 

19. 26.76 


959 


960 CANADIAN JOURNAL PHYSICS. VOL. 37, 1959 


symmetry non-zero off-diagonal matrix elements occur only between states 
having the same and therefore separated energy The samples 
have these off-diagonal elements show them small. For example 
diagonalize the truncated matrix for neutrons the eigenvalues differ 
from the diagonal values most 1.00 Mev, and the proper solution 
which belongs the eigenvalue —19.70 (cf. diagonal element 
—20.69), has The interval between the diagonal elements 
very nearly constant. Similar results hold for the other Hamiltonian matrices. 

These are admittedly only indications the actual form the self-con- 
sistent solution the H.F. equations. The effect including the complete 
set states not known. can only point out that had all the off-diagonal 
elements between occupied states been zero there would have been strong 
presumption that the original oscillator functions were self-consistent, 
that seems fair describe the results shown Tables and III nearly 
consistent. have not carried the computation any further for obvious 
reasons, but the presented results suggest that self-consistent calculation 
starting from oscillator functions would converge more rapidly than the 
similar computation for square well (Rotenberg 1955; Talman 1956). 

none these calculations does one obtain good value the binding 
energy. The average binding neutron and proton each our calcu- 
lated systems listed Table and these are all clearly too small magni- 


TABLE 
Binding energies per nucleon and the mean neutron one-body potential from H.F. calculation 
Bx, Mev Bp, Mev Mev 
—1.908 —0.583 41.00 
—3.424 —2.468 
—2.835 —2.472 


tude. cannot expect obtain correct binding energies calculation 
which ignores correlations motions nucleons nucleus (Jastrow 1955; 
Aviles 1958), and deal with those internucleon potential which describes 
interactions high energy required. However, the present computation 
was made examine the relative usefulness the oscillator wave functions 
for the initial unperturbed shell model functions, not the expectation 
properly computing the nuclear binding energy. 

the other hand, the spacing single-particle levels quite reasonable, 
and this perhaps more significant for our purpose. Tables and show 
the single-particle model value Comparison with the diagonal H.F. 
matrix elements shows that this reproduced within small error except near 


the top the system levels. 

Moreover the value derived here nearly consistent with observations 
physical nuclei. has been found (Newton 1956) that the observed level 
spacings nuclei can fitted assuming that the spacing between 
j-shells one type nucleon uniform and has, the top the occupied 
set, the value Mev. Since theré are j-shells the Ath 
oscillator shell, this yields 
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Mev. 


When the system contains closed oscillator shells the largest value 
related the neutron number (36a) that 
Thus from the level-spacing study one obtains the approximation 


(40) 
The set (39) parameters for the cutoff oscillator model yield 


The model nucleus thus has spacings very similar 
those deduced from the spacings nuclear levels. This better result than 
the calculated binding energy. However, must pointed out that equation 
(40) derived from the single-particle level density the top the occupied 
set levels while has been observed that the H.F. levels are more closely 
spaced the top. 

The H.F. study has been confined states occupied when the model nucleus 
unexcited. The degree consistency exhibited the oscillator functions 
therefore demonstrated only the region space which the ground 
state particle density appreciably different from zero. order conclude 
that the equivalent one-body potential approximately quadratic from 
the nuclear boundary require the results Section which show 
that the oscillator functions are nearly exact solutions potential which 
quadratic for and zero outside the nuclear 

obvious that this calculation leads very definite conclusion. 
have only shown indication that the oscillator single-particle wave functions 
are close the self-consistent solution for the potential (17). However, this 
statement, such is, holds for all the systems investigated. far the 
evaluation concerned, the significant measure the range the 
internucleon potential where the Yukawa constant (equation 17) 
and (u*w)/h. Sincea the “near not very sensitive 
the range the internucleon force. This was also demonstrated the 
initial calculations for Zr® various like fashion change the 
strength potential (17) will simply change proportionately. 
This will not change the regularity the interval between diagonal elements, 
though will change the relative value the off-diagonal elements. Thus the 
“near not very sensitive the strength the potential (17). 
this connection one can note that would require increase the 
strength constant bring the mean binding energies into qualitative agree- 
ment with measurement. 

The parameters (39) were derived quite arbitrary fashion. would 
clearly possible slightly increasing and retain for the model 
value close (40) and reasonable value (8.5 Mev) for the binding 
energy per neutron. This might supposed represent the results 


better variational procedure. 
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easily seen from equation (35) that the model yields qualitatively 
reasonable value for the mean binding will also reasonable value 
the separation energy the last particle, i.e. the energy required remove 
neutron 


that 


which defined (35) and the last summand the right side 
(35). The superscripts denote the nucleus which the mean binding 
nearly so. The difference (42) increases more less uniformly and slowly 
the change (42) about 2.5 Mev from oxygen lead, and this variation 
significantly increased when and are Nilsson’s (1955) values. The 
direction and order magnitude this variation agree with observation. 

thus conclude that shell model based the cutoff oscillator well 
could defined yield reasonable account the binding energies 
well exhibit the other well-known resemblances physical nuclei. 
However, not possible uniquely define the parameters the well 
these requirements. This definition should made method based 
known internucleon interactions. The partial success the present simple 
calculation allows hope that this can done. 
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APPENDIX 
illustrate the concept pseudo-potential used Section consider the 
potential consisting square attractive well with repulsive core; 


(Al) 


The corresponding pseudo potential will constructed momentum- 
dependent square well with repulsive core. For specified angular momen- 
tum between the fragments assume 


The solutions the Schrédinger equations with potentials (Al) and (A2) 
are then compared the boundary surface chosen that the 
radial logarithmic derivatives the solutions this surface have the same 
value when and the function defined that these remain 
equal each positive value One sees the momentum dependence the 
fact that the pseudo potential (A2) varies with angular momentum and 
energy. Now one can suppose that may expanded power 


and reduces the solution the transcendental equations 
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the spherical Bessel function and the combination spherical Hankel 
functions which vanishes and satisfies the Schrédinger equation con- 
taining the potential 

The depths and effective masses are listed Table 


functions energy. The case computed 184.8 Mev; 0.4 fermi 


TABLE 


E=0 87.39 V*(1) 152.6 Vo(2) 177.5 183.4 Vo(4) 184.6 


Mev M*(1)/M M*(2)/M M*(3)/M M*(4)/M 

Note: V = «~, r <0.4 fermi; V = —184.8 Mev, 0.4 fermi <r < 1.2 fermis; V = 0, r > 1.2 fermis. 
and 1.2 fermis. obvious that only the and waves are affected 


the central core and that constant effective mass for each angular momentum 
good approximation for energies 150 Mev. The and higher 
partial waves are very little the core; they are scattered mainly 
the potential step Consequently the effective mass very nearly 
equal the true mass, and the step the pseudo potential very nearly 
equal that the potential being matched. 


Pseudo potentials form (A2) with and M*/M shown match the potential 
>. 


NOTES 


ASYMPTOTIC EXPANSION OCCURRING THE 
EVALUATION LATTICE SUM* 


Recently (Maradudin and Weiss 1959) have published method for 
evaluating lattice sums and, particular, have applied sums the 
form 


a 


One the crucial steps the method applied this problem consists 
finding small expansion for the function 


n= 0 


used Poisson summation formula which required deal with 
asymptotic properties hypergeometric functions. alternative method 
which slightly simpler, and which can used even when the Fourier 
transform appearing the Poisson summation formula cannot evaluated 
exactly, use the Euler summation formula (Knopp 1951) for this problem. 
Thus find 


But can also write 

(4) 


The derivative terms can evaluated changing new variable 
which implies that (d/dx) (d/dy). This yields 


*This research was supported the United States Air Force through the Air Force Office 
Scientific Research, Air Research and Development Command under Contract No. 

tDepartment Applied Mathematics, Weizmann Institute Science, Rehovot, Israel. 

for Fluid Dynamics and Applied Mathematics, University Maryland, College 
Park, Maryland. 
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fon) 
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But known (Erdelyi 1953) that 
and 
(2k)! 
where the degree Hermite polynomial. 
Hence 


given. 

alternative, equivalent, derivation this result proceeds 
representation obtained with the aid the Poisson summation formula 
(ref. eq. (10)): 


proved Lighthill (1958) that and all its derivatives exist 


ordinary functions for and are well behaved infinity, then 


»(2m—1) 


With this result eq. (9) becomes 
at? (29) 2m dx 2m—1 (xe )e =0 


(oa) 2 1 
1 Bom d ( 


F(t) 


which just eq. (3). 
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